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ABSTRACT 

Heat-transfer  measurements  were  made  for  filmwise 
condensation  of  steam  on  externally  enhanced  horizontal 
tubes  under  vacuum  and  at  atmospheric  pressure.  Data  were 
obtained  for  copper  tubes  with  circular  fins  of  rectangular, 
triangular,  trapezoidal,  and  parabolic  cross  sections,  for 
spiral  fins  of  triangular  cross  section,  for  commercially 
available  finned  tubes  and  for  wire-wrapped  tubes.  Four 
spirally  finned  tubes  from  each  of  Cu,  Cu-Ni,  Al,  and  stain- 
less steel  and  two  tubes  with  fins  df  rectangular  cross 
section  from  each  of  Cu-Ni  and  Al  were  manufactured  and 
tested  to  investigate  the  effect  -of  thermal  conductivity. 

Among  spirally  finned  tubes,  the  optimum  fin  pitch  was 
found  to  be  1.6  mm.  The  tubes  with  a  parabolic  fin  shape 
showed  the  best  performance  with  steam- side  enhancements  of 
4.1  and  6.2  under  vacuum  and  at  atmopsheric  pressure, 
respectively.  Enhancement  ratios  as  high  as  3.5  and  2.1 
were  obtained  under  vacuum  and  at  atmospheric  pressure, 
respectively,  for  the  commercially  available  finned  tubes. 
The  heat- transfer  performance  decreased  with  decreasing  tube 
metal  thermal  conductivity. 

For  the  wire-wrapped  tubes,  an  optimum  pitch  to  wire 
diameter  ratio  of  about  5.1  was  found,  with  steam-side 
enhancements  of  about  1.9  and  2.2  under  vacuum  and  at  atmos- 
pheric pressure,  respectively.  A  recent  theoretical  anal- 
ysis of  laminar  film  of  low-surface-tension  fluids  on 
wire-wrapped  tubes  was  modified  to  include  the  condensate 
retention  of  the  tube  due  to  the  high  surface  tension  of 
water.  Agreement  between  this  modified  analysis  and  the 
experimental  data  was  favorable. 
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I.  INTRODUCTION 

A.   BACKGROUND 

It  is  well  known  that  the  power  required  to  operate  a 
naval  vessel  at  a  given  speed  is  proportional  to  its 
displacement.  Therefore,  a  major  effort  is  necessary  to 
reduce  the  displacement  in  order  to  minimize  the  required 
power.  One  of  the  largest  components  of  a  naval  vessel  is 
the  main  condenser.  In  fact,  present-day  condensers  are 
equipped  with  smooth  tubes,  and  therefore  are  large  in  size 
and  weight.  Increasing  the  performance  or  the  effective- 
ness of  the  condenser  can  reduce  the  material  and  the 
construction  cost  and  of  course  the  weight. 

The  effectiveness  of  the  condenser  is  limited  by  the 
thermal  "resistances  of  the  water  side,  the  steam  side  and 
through  the  tube  wall.  Generally,  the  thermal  resistances 
of  the  water  side  and  steam  side  are  the  most  dominant. 
Reducing  any  one  of  these  thermal  resistances  will 
contribute  to  an  improved  overall  heat-transfer  coefficient. 
Therefore,  for  a  given  heat  duty,  this  corresponds  to  a 
smaller  and  lighter  condenser.  Improved  heat- transfer 
performance  can  be  achieved  by  enhancement  of  the  water  side 
and/or  the  vapor  side.  Enhancement  on  the  water  side  is 
possible  with  turbulence  promoters,  twisted-tape  inserts, 
and  deformation  of  the  tube  to  produce  a  "roped"  scheme, 
internal  fins  or  ribs  [1] .  The  main  disadvantage  of  water- 
side enhancement  is  the  requirement  of  increased  power  for 
pumping.  Therefore,  vapor- side  enhancement  may  promise 
better  economic  advantage,  while  the  best  advantage  may  be 
achieved  by  enhancing  both  sides  based  on  a  comprehensive 
analysis.  The  enhancement  of  the  vapor-side  can  be  achieved 
by  using  low- integral  fins,  roped  tubes  or  fluted  tubes  or 
by  applying  coatings  to  promote  dropwise  condensation. 
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While  externally  finned  tubes  have  been  used  since  the 
1940s  in  order  to  enhance  the  vapor- side  coefficient  of 
tubes  used  in  refrigeration  systems,  such  tubes  have  not 
been  used  in  steam  condensers.  The  reason  for  this  appears 
to  be  the  common  belief  that  externally  finned  tubes  could 
not  enhance  steam  condensation  mainly  owing  to  the  large 
amount  of  condensate  that  floods  between  fins  in  the  lower 
portion  of  the  tube.  Since  the  surface  tension  of  water  is 
four  times  greater  than  that  of  the  refrigerants,  a  very 
significant  proportion  of  the  tube  may  trap  water  between 
fins,  which  could  result  in  poor  heat- transfer  performance. 

The  theoretical  treatment  of  the  steam  condensation 
problem  on  horizontal  finned  tubes  is  very  difficult  due  to 
the  large  number  of  controlling  parameters,  such  as  gravita- 
tional and  surface  tension  forces,  fin  spacing,  height, 
thickness  and  shape  leading  to  three-dimensional  flow  of 
condensate.  Due  to  the  complexity  of  the  problem,  any  theo- 
retical model  requires  simplifying  assumptions  which  can 
lead  to  inadequate  results.  Therefore,  a  large  pool  of 
reliable  data,  systematically  covering  all  of  the  relevant 
variables,  is  essential  in  order  to  test  simplified  theoret- 
ical models  and/or  to  arrive  at  a  satisfactory  correlation. 

This  thesis  effort  is  a  continuation  of  research  being 
conducted  at  the  Naval  Postgraduate  School  (NPS)  under  a 
grant  from  the  National  Science  Foundation.  The  basic  test 
apparatus  has  been  constructed  by  Krohn  [2]  .  Graber  [3] 
provided  the  instrumentation,  and  took  preliminary  data  as 
the  system  experienced  problems  with  non- condensing  gases 
and  partial  dropwise  condensation  on  copper  tubes.  Poole 
[4]  made  further  improvements  on  the  apparatus  especially 
for  leak  tightness.  He  operated  the  apparatus  both  under 
vacuum  and  at  atmospheric  pressure,  and  tested  a  total  of 
six  finned  tubes,  with  different  fin  spacing,  as  well  as  a 
smooth  tube.   Unfortunately,  Poole  had  problems  owing  to  the 
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occurrence  of  partial  dropwise  condensation.  Using  this 
system,  Georgiadis  [5]  was  able  to  obtain  complete  filmwise 
condensation  on  26  copper  tubes.  The  repeatability  of  data 
obtained  by  Georgiadis  proved  the  accuracy  of  the  test 
apparatus  and  associated  instrumentation  which  was  basically 
the  same  as  that  used  by  Poole  [4]  with  some  minor  modifi- 
cations. Georgiadis  tested  a  total  of  23  finned  tubes  with 
rectangular- section  fins  and  three  smooth  tubes.  He  system- 
atically varied  the  fin  spacing  (0.5,  1.0,  1.5,  2.0,  4.0  and 
9.0  mm),  fin  thickness  (0.5,  0.75,  1.0  and  1.5  mm)  and  fin 
height  (1.0  and  2.0  mm).  Table  I  shows  the  combinations  of 
fin  dimensions  used  for  these  tubes.  Based  on  both  vacuum 
and  atmospheric  runs,  Georgiadis  reported  an  optimum  fin 
spacing  of  1.5  mm  and  an  optimum  fin  thickness  of  0.75  to 
1.0  mm.  Among  the  finned  tubes  with  a  fin  height  of  1.0  mm, 
the  tube  with  a  fin  spacing  of  1.5  mm  and  fin  thickness  of 
1.0  mm  provided  the  best  heat- transfer  performance.  This 
tube  resulted  in  a  steam-side  enhancement  (i.e.,  the  ratio 
of  steam-side  coefficient  for  the  finned  tube  to  the  value 
for  the  smooth  tube  at  the  same  heat  flux)  of  about  4  and 
5.7  for  vacuum  and  atmospheric  pressure,  respectively.  He 
found  that  the  heat- transfer  performance  was  most  sensitive 
to  the  fin  spacing,  while  the  effect  of  fin  thickness  was 
relatively  small.  Further,  he  found  that  the  performance 
increased  with  increasing  fin  height.  However,  he  showed 
that  the  ratio  Eo/Ar  (i.e.,  the  enhancement  beyond  the  area 
enhancement)  decreased  with  increasing  fin  height  (for 
example,  tube  6  with  e  =  1 . 0  mm  gave  Eo/Ar  values  of  2.13 
and  3.01  for  vacuum  and  atmospheric  pressure,  respectively, 
while  tube  23  with  e  =  2 . 0  mm  gave  values  of  1.69  and  2.25). 
It  appears  that  the  surface-tension  induced  thinning  of  the 
condensate  film  diminishes  with  increasing  fin  height. 

Continuing  with  this   investigation,  Flook  [6]  tested  19 
additional  tubes  (see  Table  I  for  details).   These  tubes 
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included  two  sets  of  four  tubes  with  fin  heights  of  0.5  and 
1.5  mm,  respectively.  In  addition,  he  studied  the  effect  of 
fin  shape  using  machined  fins  of  triangular,  trapezoidal, 
and  "parabolic"  fin  shapes,  (these  tubes  had  a  fin  height  of 
1.0  mm,  a  fin  base  thickness  of  0.5  mm  and  a  fin  spacing  of 
1.5  mm  at  the  fin  root).  Flook  showed  that  the  tube  with 
parabolic  fins  (tube  38)  outperformed  the  corresponding  tube 
with  rectangular-section  fins  (tube  17)  by  10  and  15  percent 
under  vacuum  and  at  atmospheric  pressure,  respectively.  As 
also  pointed  out  by  Flook,  this  tube  did  not  have  truly 
parabolic  fins.  Like  previous  researchers  [7,8],  Flook  also 
pointed  out  that  a  fin  shape,  such  as  parabolic,  that  has  a 
continuously  decreasing  curvature  from  fin  tip  to  fin  root 
provides  considerable  thinning  of  the  condensate  film,  thus 
resulting  in  improved  heat- transfer  performance. 

Despite  considerable  achievements  made  by  Georgiadis 
and  Flook,  the  very  complicated  nature  of  the  problem  being 
studied  demands  much  more  attention.  This  includes  more 
testing  to  study  the  effect  of  fin  shape,  the  effect  of  fin 
thermal  conductivity,  performance  of  commercially  available 
tubes  and  the  enhancement  that  can  be  achieved  by  wire- 
wrapping  smooth  tubes . 

B.   OBJECTIVES 

The  main  objectives  of  this  thesis  are  as  follows: 
Take  data  on  a  number  of  tubes  to  check  the  repeat- 
ability with  previous  data  [5,6], 

Take  data  on  tubes  with  fins  of  different  shapes  (trian- 
gular, trapezoidal,  parabolic,  etc.), 

Take  data  on  commercially  available  tubes, 
Take  data  on   tubes  with  different  thermal  conductivity 
having  rectangular,  triangular,  and  spiral  fin  shapes,  and 

Take  data  on  wire-wrapped  tubes  with  different  spacing 
and  wire  diameter. 
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Develop  a  theory  to  predict  the  data  for  wire-wrapped 
tubes . 
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II.  PREVIOUS  INVESTIGATIONS  OF  FILM  CONDENSATION 
ON  EXTERNALLY  ENHANCED  HUR1ZONIAL  TUBES 

A.   GENERAL  OBSERVATIONS 

When  vapor  condenses  on  smooth  horizontal  tubes  in  a 
filmwise  mode,  the  condensate  flows  down  by  gravity  and  a 
continuous  film  always  exists  around  the  tube.  The  latent 
heat  released  by  the  vapor  will  eventually  be  absorbed  by 
the  cooling  liquid  that  flows  through  the  tube. 

The  condensate  film  resists  this  heat  flow  because  of  the 
low  conductivity  of  the  liquid.  The  resistance  increases 
as  the  film  thickness  increases.  At  the  top  of  the  tube, 
the  condensate  film  thickness  is  small  and  thereby  the 
resistance  is  low  and  it  increases  with  increasing 
distance  around  the  perimeter  of  the  tube.  Since  the 
thermal  resistance  of  the  condensate  limits  the  heat- 
transfer  performance  of  the  tube,  to  enhance  heat  transfer, 
it  is  necessary  to  reduce  condensate  film  thickness.  For 
horizontal  tubes,  thinning  of  the  condensate  may  be  achieved 
by  using  a  finned,  grooved  or  a  fluted  surface. 

In  1984,  Yau  et  al.  [9]  measured  the  enhancement 
provided  by  copper  finned  tubes  over  smooth  tubes  for  film- 
wise  condensation  of  steam.  Similar  experiments  by 
Wanniarachchi  et  al.  [10]  also  in  1984  confirmed  that  the 
observed  enhancements  were  greater  than  could  be  explained 
by  the  increased  surface  area  alone.  This  additional 
enhancement  may  be  a  result  of  the  surface-tension  forces 
which  act  to  thin  the  condensate  film.  The  effect  of 
surface  tension  was  first  described  by  Gregorig  [7]  using  a 
fluted  surface.  The  surface  tension  induced  a  large  pres- 
sure gradient  along  the  fin  surface.  This  induced  pressure 
gradient  can  be  explained  by  using  Figure  2.1. 
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Condensate 


Figure  2.1    Schematic  of  Condensate  Profile  on 

Unflooded  Fin. 

The  pressure  gradient  due  to  the  effect  of  the  surface 
tension  between  a  liquid  and  vapor  is  inversely  proportional 
to  the  radius  of  curvature  of  the  condensate  surface. 

The  pressure  of  the  condensate  at  point  A  is  higher  than 
the  vapor  pressure  because  of  the  convex  condensate  surface 
at  this  point.  The  condensate  surface  at  the  valley  is 
rather  flat.  This  nearly  infinite  radius  of  curvature  of 
the  condensate  surface  results  in  no  pressure  difference 
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induced  by  the  surface  tension  at  this  point.  Therefore, 
the  pressure  at  point  B  is  almost  the  same  as  the  vapor 
pressure.   These  pressures  are  given  by: 


A    v   rA 


P   =  P 
B    v 


(2.2) 


where 


P        =  vapor  pressure, 


Pa   Pg    =  liquid  pressure  at  points  A,  B,  and 

radius  of  curvature  of  t 
film  at  points  A,  and  B. 


a  t-q        =  radius  of  curvature  of  the  condensate 


At  point  A,  the  radius  of  curvature  is  small,  so  the  pres- 
sure at  point  A  is  higher  than  the  pressure  at  point  B  (see 
equations  (2.1)  and  (2.2)).  Since,  in  reality,  the  radius 
of  curvature  changes  along  the  condensate  surface,  between 
points  A  and  B,  the  pressure  within  the  condensate  film 
varies  along  the  height  of  the  fin.  The  overall  pressure 
difference  between  points  A  and  B  is  given  by  equation 
(2.3). 


^-k  (2'3) 


20 


where 


AP^B  =  pressure  difference  between  points  A  and  B. 

Since  the  radius  of  curvature  of  the  condensate  film  at 
point  A  is  very  small,  we  can  see  from  equation  (2.3)  that 
there  is  a  large  pressure  difference  between  points  A  and  B. 
This  pressure  difference  causes  the  condensate  to  flow  from 
point  A  to  point  B,  thinning  the  condensate  layer.  On  the 
other  hand,  the  flow  of  condensate  between  the  fins  depends 
on  the  ratio  of  surface  tension  forces  to  gravity  forces 
since  the  former  acts  to  retain  the  condensate  between  the 
fins  while  the  later  acts  to  drain  the  condensate.  As 
surface  tension  increases,  the  condensate  tends  to  flood  a 
larger  area  of  the  tube  in  which  the  condensate  layer  is 
thick  and  the  thermal  resistance  increases,  so  a  small  heat 
transfer  coefficient  results.  The  flooded  portion  of  the 
tube,  as  mentioned  in  section  A,  is  defined  by  the  retention 
angle,  (y)  (i.e.,  the  angle  from  the  bottom  of  the  tube  to 
the  highest  position  of  the  tube  where  the  interfin  space  is 
still  full  of  condensate).  The  retention  angle  depends  on 
the  fin  spacing,  surface  tension  and  gravity  forces,  and  the 
fin  shape.  Therefore,  on  the  one  hand,  using  fins  around  a 
smooth  tube  increases  the  condensing  area  and  thins  the 
condensate  film  along  the  fin  surface.  However,  these  bene- 
ficial effects  are  offset  by  the  flooding  that  occurs. 
Decreasing  the  retention  angle  increases  the  heat  transfer 
performance.  Therefore,  any  means  reducing  the  retention 
angle  is  beneficial.  As  mentioned  in  section  A,  one  way  to 
decrease  the  retention  angle  is  by  attaching  drainage  strips 
at  the  lower  part  of  the  tube. 
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B.   CONDENSATE  RETENTION 

In  1946,  the  first  measurements  of  condensate  retention 
were  made  by  Katz  et  al.  [11].  These  measurements  were  made 
under  static  conditions  (i.e.,  no  condensation  taking  place) 
using  water,  aniline,  acetone,  and  carbon  tetrachloride  on 
ten  different  tubes  with  fin  densities  from  276  to 
984  fins/m,  and  fin  heights  from  1.2  to  5.7  mm.  They  meas- 
ured the  retention  angle  by  visual  observation  and  by 
weighing  the  amount  of  retained  liquid.  Theoretical  treat- 
ment of  the  problem  using  the  measurement  of  surface  tension 
by  a  capillary  tube  and  by  the  pendant  drop  method  was  made 
to  develop  a  formula  to  predict  condensate  retention  as  a 
function  of  condensate  properties  and  the  dimensions  of  the 
tube.  Their  result  for  the  condensate  retention  angle,  vj/  is 
given  by  equation  (2.4): 


stn^ 


o 


■(4Df   -    2DQ  +   2s 
*    (Df*  -   D0z)s 
,4 


1  18Q 
980 


(2.4) 


where 


<j  =  surface  tension, 

p^  =  density  of  condensate, 

g  =  acceleration  of  gravity, 

Df  =  fin  Diameter, 

DQ  =  outside  diameter  of  tube,  and 

s  =  fin  spacing. 

It  was  shown  that  condensate  retention  depends  mainly  on  the 
ratio  of  surface  tension  to  liquid  density  and  on  the  fin 
spacing. 
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In  1981,  Rudy  and  Webb  [12]  measured  condensate  reten- 
tion angles  on  three  integral- fin  tubes  with  three  fin 
densities  (748,  1024,  1378  fins/m).  They  used  three 
different  fluids  (water,  R-ll,  and  n-pentane)  under  both 
static  and  dynamic  conditions.  Their  results  showed  that 
the  retention  angle  increases  with  increasing  surface 
tension  to  density  ratio  of  the  fluid.  They  also  showed 
that  the  difference  between  static  and  dynamic  retention 
angles  was  very  small.  For  water,  they  reported  that  a 
significant  portion  of  the  tube  surface  was  flooded. 

In  1982,  Rifert  [13]  reported  equation  (2.5)  for  the 
retention  angle  using  a  model  of  the  capillary  rise  height 
of  the  fluid  along  a  vertical  plate. 


\\>  =   -cos-1 


2  a  (P  -  Pf)~|  (2.5) 


r   2  o  (p  -  pf)~| 

L     Pf  8  D0  Ap  J 


where 


P   =  wetted  perimeter, 

Pf  =  fin  pitch,  and 

A^  =  Profile  area  of  the  fin. 

Later,  in  1983,  Rudy  and  Webb  [14]  developed  an  analyt- 
ical model  to  predict  condensate  retention.  They  used  two 
finned  sections,  one  in  tubular  form  and  the  other  by  split- 
ting the  tubular  section  and  unrolling  it  into  a  vertical 
plate.  They  found  that  the  vertical  rise  height  of  the 
condensate  was  the  same  for  these  two  cases.  Based  on  this 
observation,  they  modelled  condensate  retention  on  a  flat 
plate  to  express  the  same  on  the  finned  tube.  They  made  a 
simple  force  balance  on  the  free  body  of  condensate  and 
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Tube  Wall 


Retention  Angle  (y) 


Figure  2.2    Schematic  of  Condensate  Retention 

on  Finned  Tubes . 
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developed  an  expression  for  the  condensate  retention  angle 
as  given  by  equation  (2.6): 

t-  cos-ifl  -  2  a  (2  e  -  t)1  (2.6) 

L     Pf  §  e  s  D0  J 

where 

e  =  fin  height. 

Both  their  analytical  and  test  results  showed  that  the 
retention  angle  increases  with  increasing  surface  tension- 
to-density  ratio.  Experimental  results  involving  the  use  of 
water,  R-ll,  R-12,  ammonia,  and  n-pentane  were  predicted  to 
within  10  percent. 

Owen  et  al.  [15]  also  recognized  the  necessity  of 
including  the  effects  of  condensate  retention  in  the  heat- 
transfer  models.  The  main  simplifying  assumption  for  their 
model  was  that  the  condensate  retention  angle  was  indepen- 
dent of  condensation  rate,  so  there  is  no  difference  between 
a  static  test  and  dynamic  condensation.  Therefore,  they 
considered  only  a  static  analysis.  A  simple  force  balance 
between  surface  tension  and  gravitational  forces  resulted  in 
an  equation  for  the  condensate  retention  angle  as  shown 
below: 


t" 


i_!l  (2-7) 

Pf  g  s  Df  J, 


This  equation  is  the  same  as  equation  (2.6),  except  that 
equation  (2.7)  is  independent  of  fin  thickness  (t).  A  good 
agreement  between  this  equation  and  the  available  data  were 
reported  by  Rudy  and  Webb  [12] . 
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In  1983,  Honda  et  al.  [16]  performed  experiments  on 
finned  tubes  with  and  without  porous  drainage  plates  using 
R-113  and  methanol  as  working  fluids.  They  revealed  from  a 
photographic  study  that  the  static  and  dynamic  profiles  of 
the  retained  condensate  were  almost  the  same,  and,  by 
attaching  a  porous  drainage  plate,  they  demonstrated  a 
significant  reduction  in  the  retention  angle.  Considering  a 
force  balance  between  gravity  and  the  surface  tension  force 
acting  on  the  condensate,  they  made  a  theoretical  analysis 
to  predict  condensate  retention,  leading  to  equation  (2.8): 


r     4  o  cos6  "I  (2.8) 

\b   »  cos"  Ml ""    n  " 


where 

0  =  fin  tip  half-angle. 

They  reported  very  good  agreement  between  their  theory, 
their  own  data  and  other  available  experimental  data 
[11,12]. 

Yau  et  al.  [9]  measured  the  condensate  retention  angle 
using  water,  ethylene  glycol,  and  R-113  for  finned  tubes 
with  and  without  drainage  strips.  They  used  an  apparatus  to 
simulate  condensation  on  finned  tubes.  From  their  results, 
it  appears  that  a  drainage  strip  attached  edgewise  to  the 
bottom  of  a  tube  has  a  significant  effect  on  removing  the 
condensate,  so  liquid  retention  is  significantly  reduced. 
They  modified  equation  (2.8)  in  order  to  fit  their  experi- 
mental data,  and  developed  the  empirical  relation  given  by 
equation  (2.9): 


r         l  .66  o  cose  i 

*  =  cos~l  L1  -    Pf  g ;  nf  J  (z-?) 
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where 

0  =  fin  tip  half  angle. 

Continuing  with  their  investigation  on  condensate 
retention,  Rudy  and  Webb  [17],  in  1985,  modified  their 
previous  model  [14]  for  predicting  condensate  retention  on 
horizontal  tubes  with  fins  of  arbitrary  shape.  Experiments 
were  made  on  four  finned  tubes  with  fin  densities  from  748 
to  1378  fins/m  and  one  spine  tube  with  a  fin  density  of  1378 
fins/m.  The  fluids  tested  were  R-ll,  n-pentane,  and  water. 
In  addition,  they  tested  a  Thermoexcel-C  tube  with  fin 
density  of  1417  fins/m  and  R-ll  as  the  working  fluid.  As  in 
the  previous  models,  this  model  is  based  on  the  equal  capil- 
lary rise  height  for  a  tubular  section  and  another  section 
that  was  made  by  splitting  a  tube  section  and  unrolling  it 
into  a  vertical  flat  surface.  Equation  (2.10)  is  recom- 
mended to  predict   the  condensate  retention  angle  : 


\J/  =  cos"  *■ 


r       2  o  cp£  -  tb)    1 

L1    "   0o    Pf   g    [(tb  +   s)    e   -   Ap]J 


(2.10) 


where 

Pi  =  wetted  perimeter  of  fin  cross  section, 

t^  =  fin  base  thickness,  and 

A^  =  profile  area  of  fin  over  fin  cross  section. 

From  equation  (2.10),  the  retention  angle  increases  for 
an  increase  of  surface  tension  to  density  ratio  of  the 
liquid,  fin  density  or  for  an  decrease  of  tube  diameter. 
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For  the  case  of  a  horizontal  tube  with  rectangularly-shaped 
fins,  equation  (2.10)  reduces  to  equation  (2.7).  The  exper- 
imental deviation  from  the  predictive  value  of  equation 
(2.10)  was  ±  10  percent. 

C.   THEORETICAL  MODELS 

In  1948,  Beatty  and  Katz  [18]  performed  experiments  with 
propane,  n-butane,  n-pentane,  sulfur-dioxide,  methyl  chlo- 
ride, and  Freon-12  condensing  on  single  finned  tubes  with 
fin  densities  from  422  to  630  fins/m  to  obtain  the  vapor- 
side  heat-transfer  coefficient.  They  used  the  Nusselt  equa- 
tions for  condensation  on  a  horizontal  tube  and  on  a 
vertical  surface,  and  considered  the  finned  tube  to  be  a 
combination  of  two  parts,  a  horizontal  plain  tube  and 
vertical  fins.  Thus,  they  expressed  the  average  heat- 
transfer  coefficient  by  a  Nusselt-type  equation  based  on  an 
equivalent  diameter.  They  modified  the  customary  leading 
constant  (0.728)  found  in  the  Nusselt  equation  to  fit  their 
experimental  data  and  their  correlation  is  given  below: 


lBK 


=   0.689 


pj  pf  (pf  -  pv)  s^VTy        (2.u) 


x  -,!/»         A.    rl   -|»/»       .    .    n  Af  m    ti* 


Dej 


-j<_ri.r +i.3ii£[ii  '  (2.i2) 

Aeff  LD0J  Aeff  L*J 


where 


Ar   =    tt   D0   L 


(2.13) 
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Af  -  T   ^Df2  ~  Do2)  N  L 


Aeff  *  Ar  +  n  Af 


(2.14) 
(2.15) 


»  (Of  -  Dp2)  (2.16) 

x      ZT~D^ 

where 

Aeff  =  effective  area  of  finned  tube, 

Af  =  total  surface  area  of  finned  tube, 

Ar  =  surface  area  of  smooth  tube, 

De  =  equivalent  tube  diameter, 

^BK  =  averaSe  vapor-side  heat-transfer  coefficient, 

hfg  =  cspecif ic  enthalpy  of  vaporization, 

kf  =  thermal  conductivity  of  condensate, 

n  =  fin  efficiency, 

AT  =  vapor- side  temperature  drop, 

|i£  =  viscosity  of  condensate,  and 

pv  =  density  of  vapor. 

The  empirically  determined  leading  constant  (0.689)  in 
equation  (2.11)  is  only  5  percent  less  than  the  theoreti- 
cally derived  constant  (0.728)  using  Nusselt  theory.  But 
the  average  heat-transfer  coefficient  is  greater  than  that 
predicted  by  Nusselt  theory  for  a  smooth  tube  since  the 
equivalent  diameter  is  smaller  than  the  outside  diameter  of 
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the  tube.  They  claimed  a  maximum  error  of  +  7.2  percent  and 
-  10.8  percent  for  the  fluids  they  tested.  However,  they 
did  not  take  into  account  surface-tension  effects  for  thin- 
ning the  condensate  along  the  fin  height.  They  also 
neglected  condensate  retention  and  assumed  gravity  to  be 
entirely  responsible  for  the  flow  of  condensate.  Over  the 
decades  following  their  work,  many  researchers  have  found 
this  model  to  be  quite  adequate  for  low-surface-tension 
fluids  and  for  tubes  with  moderate  fin  densities  (i.e.,  for 
condensing  fluid-fin  density  combinations  to  yield  low 
retention  angles).  However,  as  the  fin  density  or  the 
surface  tension  increases,  the  model  tends  to  overpredict 
the  heat- transfer  coefficient  [12] . 

Some  years  later  in  1971,  analytical  and  experimental 
studies  of  condensation  on  horizontal  tubes  with  trapezoi- 
dally  shaped  fins  were  performed  by  Karkhu  and  Borovkov 
[19]  for  condensation  assuming  surface  tension  to  create  the 
dominant  force.  The  analytical  solutions  were  based  on  the 
following  assumptions:  1)  the  thin  condensate  film  repre- 
sents a  laminar  boundary  layer;  2)  surface  tension  causes  a 
pressure  gradient  along  the  fin  side;  3)  gravitational  and 
inertial  terms  in  the  equation  of  motion  of  the  film  along 
the  side  of  the  fins  were  small  compared  to  surface  tension 
terms  and  were  neglected;  4)  the  motion  of  condensate  in  the 
trough  area  is  laminar;  5)  condensate  drains  by  gravity  into 
the  trough;  6)  no  condensation  takes  place  on  the  flooded 
portion  of  the  tube;  and  7)  the  fin  temperature  is  constant 
along  the  height  of  the  fin.  Using  Nusselt's  basic  assump- 
tions and  the  differential  equation  of  condensate  motion 
(assuming  radial  flow  of  condensate  feeding  into  the 
interfin  space)  with  appropriate  boundary  conditions,  they 
were  able  to  obtain  the  thickness  of  the  condensate  film  in 
the  interfin  spacing  (equation  (2.17)).  In  order  to  calcu- 
late the  temperature  distribution  along  the  fin  height,  they 
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assumed  one-dimensional  heat  conduction.  Using  numerical 
methods  to  solve  the  resulting  differential  equations,  they 
found      expressions    for   the   heat-transfer   coefficient: 


Zb   =   1.6   I!    °-2   (1    -  0.35   H  -°-  3  m) 


h    = 


G   hfg 


AT  =   0.38   +  0.62   n+1   -   0.012   n 


where 


o    lf\    lip  k,  T      )     3/«*  Ro 
H   =  2.86 f — £-J_s_ 


sin3e   (1    +   tgO)1/4   cos1/1*    9   pf7^  e7/ 2  hf    3/l+ 


(2.17) 


FS   at  (2.18) 

where 

e  tt  D0  (2.19) 

S  2  +  cos9  2 


(2.20) 


[pf   hfg      o]     V4  kf    3/<*  e    3/2] 
1,4      [uf   b   AT]  l/H   k    [1   +   tg6]    [2   b  +  e    sin    6]  (2.21) 


2.22) 


0  cos  8  .  „    _ _  . 

We   -  r- -n c (2.23) 

b    (1   +   tan8)    e    pf 

where 

b      =    half   of    fin   tip   width, 

e      =    fin   height, 
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Fc  =  effective  condensate  surface, 

G   =  condensate  flow  rate, 

h   =  heat-transfer  coefficient, 

Z   =  dimensionless  depth  of  condensate 
between  fins , 

Z^  =  dimensionless  depth  of  condensate  at  fin  base,  and 

0   =  fin  semivertex  angle. 

Experiments  were  performed  with  four  different  finned  tubes 
to  condense  both  steam  and  R-113  with  slowly  moving  vapor 
and  when  the  Weber  number  (equation  (2.23))  is  greater  than 
10.  Using  the  measured  temperature  at  the  fin  root,  they 
found  equation  (2.21)  for  the  dimensionless  depth  of  the 
condensate  at  the  fin  base  within  ±  2  percent  of  the  experi- 
mental data.  Also,  they  solved  the  heat  conduction  equation 
over  the  fin  to  find  the  temperature  distribution  over  the 
fin  height  (equation  (2.19)).  They  found  the  vapor-side 
coefficients  to  be  50  to  100  percent  greater  than  that  for  a 
smooth  tube.  Further,  they  reported  that  their  predictions 
agreed  to  within  ±5  percent  with  their  experimental  data. 

In  1973,  Edwards  et  al.  [20]  reported  an  analytical 
model  for  condensation  on  circumferential  grooves  on  hori- 
zontal tubes  that  included  the  surface  tension  effect, 
gravity,  viscous,  capillary  pressure,  and  condensate  accel- 
eration during  the  flow  around  the  tube.  This  model  is 
based  on  the  following  assumptions:  1)  the  condensate  pres- 
sure is  uniform  over  any  cross  section;  2)  the  radius  of 
curvature  of  the  meniscus  in  the  flow  region  at  the  trough 
is  constant;  3)  the  heat  transfer  and  vapor  friction  on  the 
meniscus  are  negligible;  4)  the  draining  condensate  from  the 
fin  side   has  zero  velocity;   5)  the  grooves  have  small 
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height  compared  to  the  tube  radius;  and  6)   the  film  has  a 
contact  angle  to  the  fin  tip. 

Using  conservation  of  mass,  with  the  overall  heat- 
transfer  coefficient  as  a  function  of  local  pressure  differ- 
ence and  making  a  simple  force-momentum  balance  over  an 
element  of  condensate  film,  they  found  a  relationship  for 
the  local  heat-transfer  coefficient  as  given  by  equation 
(2.24): 


w   L  fcc  +  do 


1/2 

(2.24) 


where 


9o=  j   [z    (In  -  +  0.11593]2  cot  6g  (2.25) 


,  |"2  kf  tanM1/2  (2.26) 

-  L   8,  km   J 


me,,]1/2 

where 


k£  =  liquid  thermal  conductivity, 

1^  =  fin  thermal  conductivity, 

w   =  groove  width, 

0C  =  contact  angle,  and 

0   =  groove  half  angle. 

They  assumed  no  heat  transfer  through  the  flooded  portion  of 
the  tube.  Further,  they  assumed  only  circumferential  flow 
of  condensate,  thus  neglecting  any  flow  along  the  fin 
surface  in  the  radial  direction.  While  this  assumption 
could  result  in  poor  predictions,  they  did  not  provide  a 
comparison  of  their  theory  with  any  experimental  data. 
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In  1977,  Zozulya  et  al.  [21]  modified  their  previous 
model  [19]  to  find  expressions  for  the  rate  of  heat 
transfer.  Using  the  differential  equation  of  condensate 
motion  at  the  trough,  and  the  average  temperature  difference 
determined  by  numerical  methods,  they  arrived  at  equation 
(2.27)  which  gives  the  dimensionless  height  of  the  conden- 
sate at  the  interfin  spacing: 


Z  =  1.8  Fi 


0.32  (2.27) 


2    (W2   (kf    yf    AT)  3/t+  cos  5/**(j)  Do  /2    28) 

Fi   =   tb   tf   bfgd/l+   e^z   (pf   g)//4   (1   +  tan$)l/H 

b   =  thickness  at  top  of  fin, 

t^  =  thickness  at  fin  base,  and 

z   =  dimensionless  condensate  film  thickness 

in  the  interfin  spacing. 

They  compared  the  results  of  equation  (2.27)  and  available 
data  of  refrigerants  (R-ll,  R- 12 ,  and  R-21)  for  condensation 
on  finned  tubes  manufactured  of  copper,  brass,  and  steel 
with  rectangular  and  trapezoidal  fin-shapes.  Also,  they 
compared  experimental  data  for  condensation  of  R-113  and 
steam  on  different  finned  tubes.  Discrepancies  within  ±  15 
percent  were  reported. 

In  1979,  Webb  [22]  reported  a  procedure  for  the  design 
and  optimization  of  a  fin  surface  for  heat- transfer  perform- 
ance. Equations  (2.31),  (2.32)  and  (2.33)  were  recommended 
to  calculate  the  optimum  profile  given  by  equation  (2.29)  in 
order  to  maximize  the  heat- transfer  coefficient  given  by 
equation  (2.30): 
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SoPt   -   \53  +  hw   ) 


(2.29) 


-0.2 
hw  =   1.055  Fi    (   Fi  F2   )  (2.30) 

?l   =   K   (B    0,,,)  (2.31) 


25    o  2    L   AT  , 

F2   =  (2.32) 

n  g  Pf 


B   -  hf§   Pf  g  i  (2.33) 

Mf   kf    AT 

where 

h^   =    heat-transfer   coefficient, 

p      =    projected   area   of    convex   surface, 

S      =   value    of    s    at    0=    0m, 

0   =  angular  coordinate  measured, 

from  the  crest  of  convex  surface,  and 
P   =  P/2S. 

According  to  the  author,  this  model  underpredicts  the  heat- 
transfer  coefficient.  The  calculated  augmentation  ratio 
based  on  the  projected  surface  area  h_/hj,ju  ranged  from  3.4 
to  3.8  for  tubes  with  length  from  4  ft  to  40  ft  while  his 
experiment  showed  values  in  the  range  from  4  to  8. 

In  1980,  Rifert  [23]  studied  condensation  of  stationary 
vapor  on  horizontal  finned  tubes  enhanced  by  the  effect  of 
surface-tension  forces  that  tend  to  pull  the  condensate  to 
the  fin  root.   In  his  analysis,  he  divided  the  tube  into  two 
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zones:  a)  the  unflooded  zone  where  the  condensate  film  is 
thin,  and  b)  the  flooded  zone  where  the  condensate  film  is 
thick.  He  solved  the  two-dimensional  heat-conduction 
problem  for  the  wall  by  numerical  methods  for  each  zone  and 
then  determined  the  mean  heat  flux.  In  cases  where  conden- 
sate is  retained  in  more  than  half  of  the  tube  perimeter, 
Rifert  pointed  out  that  a  three-dimensional  form  of  the 
heat-conduction  equation  must  be  used.  Solutions  to  these 
equations  by  numerical  methods  revealed  that,  in  most  cases, 
the  fin  temperature  is  very  nonuniform  and  it  depends  on  the 
properties  of  the  wall  and  the  vapor  and  the  heat  flux.  He 
stated  that  for  the  highly  non-isothermal  fin  surface,  the 
use  of  and  average  temperature  drop  from  vapor  to  the  outer 
wall  temperature  (AT)  yields  computed  heat  flux  values  that 
are  very  sensitive  to  AT.  Since  this  is  unacceptable,  he 
recommended  the  use  of  the  average  heat  flux  for  the 
computation. 

In  1981,  Adamek  [24]  presented  a  method  for  the  design 
of  an  optimum  surface  for  condensing  heat- transfer  perform- 
ance. Similar  to  other  researchers  [7,8,22],  Adamek  recog- 
nized the  importance  of  surface  tension  on  the  heat- transfer 
performance  of  finned  surfaces .  Since  the  dominant  force  on 
the  crest  is  the  surface  tension,  he  neglected  gravitational 
forces  in  this  region.  He  derived  equations  for  the  conden- 
sate film  thickness  (equation  (2.34))  and  the  wall  surface 
profile  by  defining  the  curvature  as  a  function  of  the 
distance  along  the  surface  (equation  (2.35)).  Using  equa- 
tion (2.36)  for  the  the  curvature  of  the  profiles  of  the 
wave  crest,  and  the  necessity  that  the  pressure  within  the 
condensate  must  decrease  from  wave  crest  to  the  trough,  he 
defined  a  family  of  condensate  surface  profiles,  whose 
curvatures  are  given  by  equation  (2.38).  A  number  of  £ 
values  and  their  corresponding  condensate  surface  profiles 
are  shown  in  Figure  2.3.   He  found  equations  (2.37),  (2.39), 
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and  (2.40)  for  the  film  thickness,  the  average  heat- transfer 
coefficient  and  condensate  flowrate,  respectively.  The 
parameter  |  in  equation  (2.37)  characterizes  the  aspect 
ratio  of  the  fin  cross  section  (ratio  of  the  height  to  the 
thickness).  As  the  aspect  ratio  increases,  the  parameter  £ 
decreases.  As  shown  by  Adamek,  £  =  -  0.5  represents  the 
optimum  surface  for  maximum  values  of  the  condensate  flow- 
rate  and  the  average  heat- transfer  coefficient. 

6 
'  6(s)  =  |(<  r1/3  (A   /  C  (<  )i/3  ds  +  C0  1  i/*  (2.34) 


W(s)  =  f(s)  -  5(s)  nf(s)  (2.35) 


<(s)  =  a  s   -  <0       -1  <  C  <  0  (2.36) 


e+i  „  2-n   v* 


6(s)  -   12pf  Uf  M  1   S,,^  s  2-H  (2i37) 


l   em  (c  +  l)  f  _  /s_nC"| 
r   Sm    5    L    ^sm/  J 


(2.38) 


h  -  2  149  kf  r°  hf^  gs  Sm  Pf  U  +  'M17'  (2-39) 

sn  L     Mf  kf  at        TTTTpJ 
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(2.40) 


where 

W(s)  =  wall  profile, 

f(s)  =  film  profile, 

k(s)  =  local  curvature  of  the  condensate  surface, 

s  =  length  of  path  in  liquid  film, 

S     =  length  of  convex  surface  over  which 
the  condensate  flow, 

8  =  film  thickness, 

0  .=  rotation  angle  of  normal  to  fin  surface, 

0m  =  angle  from  origin  to  Sm, 

£  =  ratio  of  slenderness,  and 

r  =  radius  of  curvature. 

In  1981,  Shklover  et  al.  [25]  treated  film  condensation 
for  finned  tubes  to  investigate  the  effect  of  metal  thermal 
conductivity  on  the  heat- transfer  performance.  Stationary 
steam  was  used  as  the  condensing  fluid.  They  showed  that 
as  the  thermal  conductivity  decreases,  the  temperature 
difference  through  the  film  decreases  and  the  temperature 
difference  through  the  wall  increases.  For  this  reason,  the 
finned  tubes  made  of  stainless  steel  or  german  silver  have 
the  same  heat- transfer  performance  as  the  smooth  tube,  while 
brass  and  copper  finned  tubes  outperformed  the  smooth  tubes. 

Rudy  and  Webb  [12]  proposed  a  possible  improvement  to 
the  Beatty  and  Katz  [18]  model  by  taking  into  consideration 
condensate  retention.   They  applied  equation  (2.11)  only  for 
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Figure  2.3    Adamek  [24]  Condensate  Surface  Profiles. 

the  unflooded  part  of  the  tube  and  recommended  equation 
(2.41)  for  the  average  heat- transfer  coefficient. 


h  =  hBr 


^1 


(2.41) 
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where  hfiK  is  computed  using  equation  (2.11).  But,  since 
this  is  a  gravity-based  model  and  it  neglects  any  heat 
transfer  through  the  flooded  portion  of  the  tube,  it  under- 
predicts  the  average  heat-transfer  coefficient  of  condensing 
R-ll  by  10  to  30  percent.  They  recognized  that  the  surface 
tension  effect  must  be  taken  into  account. 

In  1982,  Webb  et  al.  [26]  developed  a  new  model  which 
included  surface  tension  effects.  They  modified  the  original 
Nusselt  equation  for  a  vertical  plate  so  that  surface 
tension  causes  the  condensate  to  drain  from  the  fin  tip  to 
the  base  and  gravity  causes  the  condensate  to  flow  in  the 
channel  between  the  fins.  Assuming  surface  tension  as  the 
dominant  force  along  the  fin  side,  they  proposed  equation 
(2.42)  for  the  fin  side  coefficient  hf^n  and  the  average 
heat-transfer  coefficient  for  the  entire  tube  is  given  by 
equation  (2.43): 


hfin-  o'.9«[%^f[£(i+if         (2.42) 

where  h^  was  computed  using  the  Nusselt  [27]  equation,  and 


„  =  -^hb  +  nf  -^h£ln  (2.43) 

This  model  predicted  the  heat-transfer  coefficient  within 
±  10  percent  for  R-12. 

Using  equation  (2.7)  for  the  condensate  retention  angle 
they  developed  earlier,  Owen  et  al.  [15]  modified  the 
Beatty  and  Katz  model  to  include  the  retention  angle.  They 
divided  the  tube  into  two  parts:  the  unflooded  portion,  and 
the  flooded  portion  with  the  condensation  occuring  on  both 
the  retained  condensate  and  the  fin  tips.  The  equations 
necessary  for  this  model  are  listed  below: 
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h-(,~  gh.+  *i,f  £■  (2.44) 

where  1^  was  computed  using  equation  (2.11)  with  a  leading 
coefficient  of  0.725  instead  of  the  value  of  0.689,  and 

1      1   -1 
hf  "  (  h—  +  7-   )  2.45) 

neff    nc 

where  hc  was  computed  using  Nusselt  [27]  equation,  and 


leff 


'111  (2.46) 


keff  "  (1  "  sN)  kfin  +  sNkf  (2.47) 


where 


1^    =  heat-transfer  coefficient  for  the 
unflooded  portion  of  the  tube, 

h    =  heat- transfer  coefficient  for  a 
plain  tube, 

hf        =    heat-transfer  coefficient  for  the 
flooded  portion  of  the  tube, 

h  ££  =  heat-transfer  coefficient  of  the 

combined  fin  and  retained  condensate, 

k  ff  =  effective  thermal  conductivity, 

kf^n  =  fin  thermal  conductivity, 

N  =  number  of  fins  per  meter,  and 

\j/  =  retention  angle. 
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s    =  fin  spacing 

Owen  et  al.  showed  their  model  to  predict  the  data  for  R-ll, 
R-22,  methyl  chloride,  n-pentane,  sulfur  dioxide,  propane, 
and  n-butane  to  within  ±  30  percent.  However,  as  shown  by 
Honda  and  Nozu  [16] ,  this  model  overpredicted  the  steam  data 
by  up  to  a  factor  of  2. 

In  1983,  Rudy  and  Webb  [29]  developed  a  model  based  on 
surface-tension-induced  linear  pressure  gradient  along  the 
fin  height,  thus  assuming  radial  flow  of  the  condensate 
feeding  into  the  interfin  space.  Further,  they  assumed 
gravity-drained  flow  of  condensate  in  the  space  between 
fins.  The  Nusselt  equation  for  horizontal  tubes  was  used 
for  the  tube  area  between  fins,  while  the  fin  surface  was 
treated  by  replacing  the  body-force  term  (i.e.,  "pg")  in  the 
Nusselt  equation  by  an  equivalent  expression  based  on 
surf ace- tension- induced  pressure  gradient  as  developed  by 
Webb  et  al.  [26]  and  Rudy  [28]  earlier.  Once  again,  they 
assumed  no* heat  transfer  through  the  flooded  portions  and 
the  resulting  expression  is  given  by  equation  (2.48): 


h  = 


0.725 


+   0.943  n 


(2.48) 


i/t 


(*  -  *) 


where 


L    =  length  of  tube, 

A^t  =  surface  area  of  tube  between  fins,  and 


'ft 


=  fin  surface  area. 
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This  expression  provided  an  accuracy  of  better  than  10%  for 
condensation  of  R-ll  on  short,  finely-spaced  fins,  but  the 
accuracy  dropped  sharply  with  fins  of  increasing  height  and 
for  larger  fin  spacing.  This  was,  according  to  the  authors, 
due  to  the  assumed  linear  pressure  gradient  on  the  fin 
surface  as  this  model  is  not  valid  when  gravity  forces 
become  dominant  (i.e.,  as  e  increases).  Therefore,  equation 
(2.48)  is  valid  for  fin  densities  from  1200  to  1400  fins/m, 
and  fin  heights  of  less  than  1  mm. 

Continuing  their  research  on  film  condensation,  Honda  et 
al.  [16]  did  experiments  on  horizontal  finned  tubes  by 
attaching  a  vertical  drainage  strip  at  the  bottom  of  the 
tube  to  reduce  condensate  retention.  Using  R-113  and  meth- 
anol as  condensing  fluids,  they  found  vapor-side  enhancement 
ratios  (compared  to  the  case  without  drainage  strips)  as 
high  as  1.36  for  R-113  and  2.08  for  methanol. 

In  1984,  Honda  and  Nozu  [30]  developed  an  analytical 
model  for  film  condensation  on  horizontal  low  integral-fin 
tubes.  They  divided  the  tube  into  flooded  and  unflooded 
regions.  This  model  is  based  on  the  following  assumptions: 
1)  the  wall  temperature  is  uniform;  2)  the  flow  is  laminar; 
3)  the  condensate  film  thickness  is  small;  4)  the  dominant 
flow  on  the  fin  is  in  the  radial  direction.  Based  on  these 
assumptions,  expressions  for  Nusselt  number  representing  the 
flooded  and  unflooded  regions  were  found.  The  average 
Nusselt  number  is  given  by  equation  (2.49): 

Nud  =  (Nudu  nu(l  -  TW)H   +  Nudf  nf(l  -  Twf)(l  -  *f)       ^  ^ 
/(I  "  W^f  +  (1  -  Twf)(l-*f) 

where 

n^j    =  Fin  Efficiency, 
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ri£    =  fin  efficiency  for  the  flooded  region, 

Nuj   =  average  Nusselt  number, 

Nuj   =  average  Nusselt  number  for  the  unflooded  region, 

Nujc  =  average  Nusselt  number  for  the  flooded  region, 

q>    =  angular  coordinate, 

Twu   =  dimensionless  temperature  for  the  unflooded 
region,  and 

Tw£   =  dimensionless  temperature  for  the  flooded  region. 

Comparison  of  the  results  of  this  model  with  the  available 
experimental  data  showed  agreement  to  within  ±  20  percent 
for  11  fluids  and  22  finned  tubes.  However,  their  model 
overpredicted  steam  data  by  up  to  40  percent. 

In  1985,  Rudy  and  Webb  [31]  modified  their  previous 
models  taking  into  account  surface-tension  effects  on  film 
drainage  and  condensate  retention.  They  treated  the  conden- 
sation problem  considering  two  major  regions:  unflooded  and 
flooded  regions.  They  further  divided  the  unflooded  region 
into  finned  area  and  the  interfin  area.  They  computed  the 
average  heat- transfer  coefficient  for  the  entire  tube  as 
given  by  equation  (2.29). 


hRW 


h  n0  -  [b,  £  ♦  .  hfln  £]fLlJ}-  hf  i  (2.50) 


In  order  to  compute  the  heat- transfer  coefficient  for  the 
finned  area  (nf in) >  they  used  an  expression  developed  by 
Adamek  (equation  2.39)).    One  of  the  profiles  they  used  for 
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the  trapezoidal- shaped  fin  is  shown  in  Figure  2.4.  They 
used  equation  (2.38)  with  0m  =  85  degrees  and  for  each  fin 
geometry  an  iterative  procedure  to  establish  the  P  value  for 
each  profile  to  correct  the  fin  thickness  at  the  fin  base 
given  by  equation  (2.51). 


t  =  tfc  -  tt  +  2  5(Sm)  (2.51) 

Further,  since  they  assumed  that  the  length  of  the  convex 
surface  is  from  fin  tip  to  fin  base,  they  corrected  for  the 
film  thickness  equation  (2.52)  resulting  from  the  additional 
condensation  at  the  fin  tip: 


. "Sfo-  %+  t/2  -  5r  (2-52) 

where 


H*L  (2.53) 

where  hr  is  the  heat-transfer  coefficient  for  the  interfin 
area  in  the  flooded  region.  In  order  to  calculate  the 
interfin  area  (hr),  they  used  the  Nusselt  equation  with  an 
iteration  to  account  for  the  additional  condensate  drainage 
from  the  fins.  Finally,  to  compute  the  heat- transfer  coef- 
ficient in  the  flooded  region  they  used  the  following 
equations : 

6  =  qh2  /  Ibl  (2.54) 


(2.55) 

AT 


,     =  k  _   -  hf  AT 
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=  6-1 

e 


(2.56) 


where  qui  is  the  heat  flux  if  the  fin  thickness  were  zero. 
Using  a  numerical  technique,  they  computed  q^  based  on 
two-dimensional  conduction  through  the  fins  and  the  conden- 
sate film. 
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Figure  2.4    Fin  Geometry  for  the  Webb  et  al.  Model  [31] 
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D.   FILMWISE  CONDENSATION  ON  WIRE- WRAPPED  TUBES 

Similar  to  fin  tubes,  surface-tension  effects  can  be 
beneficial  for  filmwise  condensation  on  wire-wrapped  tubes. 
However,  the  surface-tension  effect  on  wire-wrapped  tubes  is 
different  than  that  on  finned  tubes  mainly  because  of  the 
very  little  heat  transfer  through  the  wires  compared  to  fins 
that  would  transfer  the  majority  of  the  heat. 

As  can  be  seen  from  Figure  2.5,  the  condensate  surface 
in  the  space  between  wires  on  the  wire-wrapped  tube  is 
rather  flat  and  the  pressure  difference  between  the  conden- 
sate and  the  surrounding  vapor  is  zero.  However,  due  to 
the  existence  of  a  concave  condensate  surface  at  the  point 
of  contact  between  the  wire  and  the  condensate  surface,  a 
reasonable  pressure  difference  within  the  condensate  will 
exist  from  the  inter-wire  space  to  the  immediate  vicinity  of 
the  wire.  In  fact,  the  pressure  at  point  A  is  the  smallest, 
thus  resulting  in  a  condensate  flow  toward  the  wire,  at 
which  point  the  condensate  would  rapidly  flow  around  the 
tube  and  eventually  leave  the  tube  at  the  lowest  point. 
This  axial  flow  of  condensate  toward  the  wires  generally 
results  in  a  smaller  film  thickness  than  on  a  plain  tube. 
On  the  other  hand,  the  film  is  generally  quite  thick  in  the 
region  between  the  wire  and  the  tube,  resulting  in  negli- 
gible heat  transfer  through  this  region.  As  can  be  seen  from 
Figure  2.5,  the  radius  of  curvature  of  the  condensate 
surface  near  the  wire  strictly  depends  on  the  wire  diameter. 
Further,  the  extent  of  thinning  of  the  condensate  film 
depends  on  the  wire  spacing.  Based  on  the  above-mentioned 
observations,  an  optimum  combination  of  wire  diameter  and 
pitch  must  exist  to  yield  the  best  heat- transfer 
performance . 

In  1985,  analytical  and  experimental  studies  of  conden- 
sation on  horizontal  wire-wrapped  tubes  were  performed  by 
Fujii  et  al.  [32].  Their  model  is  based  on  the  following 
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assumptions:  1)  the  liquid  film  is  very  thin  compared  to 
the  diameter  of  the  tube;  2)  the  condensate  film  flow  is 
laminar;  3)  the  inertia  forces  in  the  momentum  equation  of 
the  condensate  are  negligible;  4)  the  shear  stress  at  the 
vapor-liquid  interface  is  negligible;  5)  the  properties  of 
the  condensate  are  constant;  and  6)  the  temperature  of  the 
cooling  surface  is  uniform.  They  also  assumed  that  heat 
transfer  occurred  only  through  the  thin  film  between  wires 
and  they  neglected  any  heat  transfer  through  the  wires. 
They  developed  semi- theoretical  equations  to  predict  a  heat- 
transfer  coefficient  enhancement  ratio  as  shown  below: 


Niy  _    s     F2(A) 
Nus   (s  +  Dw)   F2(0) 


(2.57) 


where 


F2(A)  -  ^  or   Fi(<j>,A)  <H 


(2.58) 


and 


(tan(<})/2))A/3  (sinc^)1/3 
JT(tan(<{>/2))     (sin$)    d  <p 


'] 


(2.59) 


where 


A  = 


A  o  D, 


Pf  g  s^  rs 


(2.60) 


(2.61) 
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They  noticed  that  for  A  >  15 


F^o.A)  =  Fi(0,A)  =  ( 


4  (1  +  A)  \  i/" 


Therefore,  in  this  situation  with  A  >  15, 


(2.62) 


p»  (1  +  A)l  l/.*» 
*2<A)  =[ 3 J 


ih 


^2(0)  *   (3  ) 
and  equation  (2.57)  simplifies  to1  : 


Nus    (s  +  Dw)[_    3    J 


x- 


? 


(2.63) 


(2.64) 


(2.65) 


They  selected  a  value  of  0.03  for  the  coefficient  C  in  equa- 
tion (2.61)  in  order  to  fit  their  experiment el  data. 
Experiments  were  performed  with  three  wire-wrapped  tubes  to 
condense  ethanol  and  R-ll  vapor.  Enhancement  ratios  of  3.7 
and  3.3  for  R-ll  and  ethanol,  respectively,  were  reported 
(see  Figure  2.6).  The  maximum  enhancement  ratio  occurred  at 
a  p/Dw  value  of  2. 


*In   the   original   Fuiii   et   al.    paper   equation 
(2.65)  contained  an  error  by  retaining  the  term  (4/3) w 
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Figure  2.5    Condensate  Film  Profile  on  Wire-Wrapped  Tubes 
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Figure  2.6    Comparison  Between  Experimental  Data 
and  Semi-Theoretical  Model  of  Fujii  et  al.  [32]. 
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III.  DESCRIPTION  OF  TEST  APPARATUS 

A.   TEST  APPARATUS 

The  test  apparatus  used  for  this  investigation  was 
essentialy  the  same  as  used  by  Georgiadis  and  Flook  [5,6]. 
A  schematic  of  this  apparatus  is  shown  in  Figure  3.1.  Steam 
was  generated  using  distilled  water  in  a  304.8  mm  (12  in.) 
Pyrex  glass  section  which  was  fitted  with  ten  4000-Watt, 
440-Volt  Watlow  immersion  heaters.  The  steam  from  the 
boiler  flowed  upward  and  passed  through  a  304.8  mm  (12  in.) 
to  152.4  mm  (6  in.)  reducing  section  to  a  2.44  m  (8  ft.) 
long  section  of  Pyrex  glass  piping.  The  steam  flowed 
through  a  180-degree  bend  and  entered  a  1.52-m-long  section 
before  finally  entering  the  stainless-steel  test  section, 
which  is  shown  in  Figure  3.2.  The  test  tube  was  mounted 
horizontally  in  the  test  section.  A  portion  of  the  steam 
condensed  on  the  test  tube,  while  the  excess  steam  travelled 
downward  and  condensed  in  the  auxiliary  condenser.  The 
condensate  drained  back  to  the  boiler  by  gravity,  completing 
the  closed-loop  operation  of  the  system. 

The  exit  side  of  the  test  tube  was  provided  with  a 
mixing  chamber  for  accurate  measurement  of  the  outlet 
temperature  of  the  coolant.  A  view  port  was  provided  in 
the  test  section  to  allow  visual  observation  of  the  conden- 
sation mode  to  ensure  complete  filmwise  condensation  during 
data  collection.  The  auxiliary  condenser  consisted  of  two 
9.5  mm  (3/8  in.)  diameter  water-cooled  copper  tubes  heli- 
cally coiled  to  a  height  of  457  mm  (18  in.).  The  auxiliary 
condenser  was  cooled  by  a  continuous  supply  of  tap  water 
through  a  flow  meter.  A  throttle  valve  was  provided  to 
control  the  flow  rate  through  the  auxiliary  condenser,  thus 
keeping  the  system  at  the  desired  internal  pressure.  For 
example,   when   the   flowrate   through   the   test   tube  was 
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decreased,  the  flowrate  through  the  auxiliary  condenser  had 
to  be  increased.  Filtered  tap  water  was  collected  in  a 
large  sump  with  a  capacity  of  about  0.4  cubic  meters  (Figure 
3.3),  and  was  used  to  cool  the  test  tube.  Two  centrifugal 
pumps  connected  in  series  took  the  water  from  the  sump  and 
passed  it  through  a  flow  meter  into  the  test  tube.  A  valve 
on  the  discharge  side  of  the  second  pump,  and  before  the 
flow  meter,  allowed  the  velocity  of  water  flowing  through 
the  test  tube  to  be  varied  from  0  to  4.4  m/s  (14.4  ft/sec). 

A  vacuum  pump  was  operated  continuously  during  the  oper- 
ation of  the  apparatus  to  remove  non- condensing  gases  from 
the  test  section.  The  system  used  to  remove  non- condensing 
gases  is  shown  in  Figure  3.3.  It  was  unavoidable  that  the 
vacuum  pump  mainly  drew  steam  with  trace  amounts  of  air 
(non- condensing  gases).  To  minimize  the  contamination  of 
the  pump  by  the  steam,  another  condenser  was  provided  to 
condense  as  much  steam  as  possible.  This  condenser  was 
cooled  with  the  filtered  tap  water  before  it  entered  the 
large  sump.  The  condensate  from  this  steam  collected  in  a 
Plexiglas  cylinder  to  be  drained  later. 

B.   INSTRUMENTATION 

The  electrical  power  input  to  the  boiler  immersion 
heaters  was  controlled  by  a  panel-mounted  potentiometer.  In 
order  to  compute  the  input  power  to  the  boiler,  a  root-mean 
converter  with  an  input  voltage  of  440  VAC  generated  a 
signal  which  was  fed  to  the  data  acquisition  system.  A 
more-detailed  description  of  the  boiler  power  supply  is 
provided  by  Poole  [4] .  The  temperatures  of  the  steam, 
condensate  and  the  ambient  surroundings  were  measured  using 
calibrated  copper- constantan  thermocouples  made  of 
0 . 25-mm-diameter  wires.  Two  of  them  were  used  for  the  steam 
temperature,  one  for  the  condensate  return  and  one  for  the 
ambient  temperature.  These  thermocouples  had  an  accuracy 
within  ±  0.1  K.   when  compared  against  a  platinum-resistance 
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thermometer.  Since  the  temperature  rise  of  the  coolant 
through  the  test  tube  is  the  most  critical  measurement, 
considerable  attention  was  paid  to  obtaining  the  highest 
possible  accuracy.  For  this  purpose  two  independent  temper- 
ature measurement  techniques  were  used:  a  Hewlett-Packard 
(HP)  2804A  quartz  thermometer  with  two  probes  having  an 
accuracy  of  ±  0.02  K,  along  with  a  10-junction,  series- 
connected  copper- constantan  thermopile  with  a  resolution  of 
0.03  K. 

For  most  of  the  data  collected,  the  quartz  thermometer 
and  the  thermopile  agreed  to  within  ±  0.03  K  and  when  the 
difference  was  more  than  ±  0.05  K,  the  data  set  was  disre- 
garded and  a  repeat  set  was  made.  The  cooling  water  flow 
rate  was  measured  using  a  calibrated  rotameter  and  the  value 
was  fed  manually  to  the  computer.  Another  rotameter  was 
provided  to  measure  the  cooling  water  flow  rate  through  the 
auxiliary  condenser. 

A  pressure  tap  located  about  50  mm  above  the  test  tube 
was  connected  to  a  U-tube,  mercury-in-glass  manometer  gradu- 
ated in  millimeters  to  measure  the  absolute  pressure  of  the 
system.  At  the  beginning  and  at  the  end  of  each  test  run, 
an  accurate  pressure  reading  was  made  and  entered  into  the 
computer.  The  measured  system  pressure  and  the  saturation 
pressure  corresponding  to  the  measured  steam  temperature 
were  used  to  compute  the  concentration  of  any  air  that  might 
have  been  present.  For  this  purpose,  a  Gibbs-Dalton-type 
relationship  was  used.  The  computed  non- condensing  gas 
concentration  was  found  to  be  within  -  1.5  to  0  percent. 
Such  a  value  revealed  that  major  air  leaks  did  not  take 
place  following  the  last  vacuum  test  on  the  apparatus. 

C.   VACUUM  INTEGRITY 

Vacuum  tightness  for  any  condensation  system,  especially 
at  low  pressures  similar  to  marine-vechicle  condensers  which 
operate  at  an  absolute  pressure  of  2  inHg,  is  very 
important . 
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The  reason  for  this  is  because  any  small  amount  of  air  or 
other  non- condensing  gas  present  with  the  condensing  vapor 
tends  to  accumulate  at  the  liquid-vapor  interface.  When 
this  phenomenon  takes  place,  an  added  thermal  resistance 
occurs  at  the  interface,  which  will  degrade  the  heat- 
transfer  performance  considerably.  Therefore,  in  order  to  be 
able  to  collect  consistent  and  reliable  data,  exteme  care 
was  taken  to  ensure  a  leak-tight  apparatus.  In  fact,  during 
the  early  stages  of  this  investigation,  a  major  leak  was 
found  through  the  screws  at  the  test  section.  After  this 
was  fixed,  a  vacuum  test  was  carried  out  routinely  about 
once  a  month.  A  leak  rate  that  corresponds  to  a  pressure 
rise  of  about  2  mmHg  in  24  hours  nearly  at  the  operating 
pressure  was  found.  Also,  since  the  vacuum  pump  was  oper- 
ated continuously  during  the  experiment,  any  accumulation  of 
non- condensing  gases  within  the  apparatus  was  effectively 
eliminated. 

D.  DATA  ACQUISITION  SYSTEM 

An  HP-9826A  computer  was  used  to  control  an  HP  3497A 
Data  Acquisition  System  to  monitor  the  system  temperatures 
and  boiler  input  power  (using  the  converter  signal).  Raw 
data  were  processed  immediately  and  stored  on  diskette  for 
reprocessing  at  a  later  time.  After  all  the  runs  were 
collected,  the  data  were  reprocessed  using  a  new  Sieder-Tate 
Coefficient  found  by  the  modified  Wilson  method. 

E.  TUBES  TESTED 

For  this  thesis  effort,  a  total  of  twenty  six  tubes  were 
manufactured.  Table  I  lists  all  the  finned  tubes  tested  and 
their  dimensions,  including  two  tubes  tested  also  by 
Georgiadis  [5]  and  Flook  [6]  and  four  tubes  tested  by  Flook. 
Twenty  one  of  them  were  made  from  copper,  two  from  copper- 
nickel,  two  from  aluminum  and  one  from  stainless  steel. 
Each  tube  had  a  133.4  mm  length  exposed  to  steam.   Tubes  45, 
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46,  47,  52  thru  56  were  finned  tubes  with  bore  diameter  of 
12.7  mm  and  fin  root  diameter  of  a  19.05  mm.  Tubes  49,  51, 
and  57  thru  62  were  finned  tubes  with  fin  root  diameter  of 
13.7  mm.  Table  II  lists  all  the  wire-wrapped  tubes  tested 
and  their  dimension.  Tubes  63  thru  71  had  smooth  exteriors 
with  a  bore  diameter  of  12.7  mm  and  an  outside  diameter  of  a 
19.05  mm.  Tube  36  and  tubes  45  through  47  consisted  of  a 
family  of  spiral  tubes,  with  triangular- shaped  fins. 
Figure  3.8  shows  a  photograph  of  these  four  tubes  which  had 
a  fin  height  of  1.0  mm  and  pitch  of  1.06,  1.6,  2.1,  and  2.5 
mm  (tubes  36,  45,  46,  and  47). 

Two  commercially  available  tubes,  manufactured  by  High 
Performance  Tubes,  Inc.,  were  tested  to  investigate  their 
heat- transfer  performance  compared  to  the  other  machined 
finned  tubes.  They  were  finned  tubes  (tubes  49  and  51)  with 
fin  height  of  a  0.75  mm  and  fin  density  of  1102  and  1181 
fins/m,  respectively.  A  smooth  tube  (tube  50)  was  also 
prepared  by  machining  off  the  fins  so  that  the  effect  of 
fins  can  be  determined.  These  three  tubes  had  the  same 
outside  root  diameter  (17.5  mm).  In  addition  to  these 
tubes,  to  study  the  effect  of  fin  shape  on  the  heat- transfer 
performance,  four  finned  tubes  were  manufactured  to  complete 
two  sets  of  tubes  with  different  fin  profiles.  The  first 
set  consists  of  the  tubes  06,  54,  55,  and  56  with  rectan- 
gular, parabolic,  trapezoidal,  and  triangular  fin-shapes, 
respectively,  with  a  fin-base  thickness  of  1.0  mm.  The 
other  set  consists  of  the  tubes  17,  38,  53,  and  52  with 
rectangular,  parabolic,  trapezoidal,  and  triangular  fin- 
shapes,  respectively,  with  a  fin-base  thickness  of  0.5  mm. 
Each  tube  in  these  two  sets  has  the  same  fin  height  of 
1.0  mm  and  fin  base  spacing  of  1.5  mm.  Figures  3.4  shows  a 
photograph  of  tubes  06,  54,  55,  56,  while  Figure  3.5  show 
schematic  cross-sectional  of  the  "parabolic"  fin.  Figures 
3.6  and  3.7  show   a  photograph  and  cross-sectional  views  of 
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tubes  17,  38,  52,  and  53,  respectively.  As  can  be  seen  from 
Figure  3.5  and  3.7  these  tubes  do  not  have  the  exact  shapes 
as  stated  above.  For  example,  careful  examination  of  the 
schematic  cross-sectional  of  the  parabolic  tube  (tube  54) 
showed  that  the  fin- shape  is  almost  straight  near  the  fin 
base  and  circular  at  the  fin  tip,  while  the  parabolic  fin  in 
Figure  3.7  shows  that  it  had  almost  straight  sides  near  the 
fin  base  with  a  sharp  leading  edge.  Also,  as  seen  in  the 
same  figure,  a  distinction  between  triangular  and  trape- 
zoidal fins  is  not  possible.  The  reason  for  this  is  that 
the  very  thin  fins  lead  to  nearly  the  same  fin  thickness  at 
their  tips,  because  of  the  difficulties  associated  with  the 
machining  process. 

In  order  to  test  the  effect  of  fin-metal  thermal  conduc- 
tivity on  the  heat- transfer  performance,  four  spiral  tubes 
with  triangular  fin  profile  were  manufactured:  one  each  from 
copper,  copper-nickel  ,  stainless  steel  and  aluminum  (tubes 
57,  58,  59  and  60,  respectively).  In  addition,  two  tubes 
with  rectangular  fin  profiles  from  copper-nickel  and 
aluminum  were  manufactured  (Tubes  61  and  62,  respectively). 
Due  to  the  low  thermal  conductivity  of  these  tubes,  a 
smaller  outside  diameter  (13.5  mm)  was  selected  to  minimize 
the  tube  metal  resistance  (note  that  these  tubes  have  the 
same  nominal  inside  diameter  as  the  other  tubes).  The 
spiral  tubes  had  a  fin  base  thickness  of  2.1  mm,  while  the 
rectangular  fins  had  a  fin  thickness  of  1.0  mm  and  a  fin 
spacing  of  1.5  mm.  Figure  3.8  shows  a  photograph  of  these 
tubes . 

Finally,  nine  smooth  copper  tubes  were  manufactured,  and 
each  was  wrapped  with  a  titanium  wire  (0.5,  1.0  or  1.6  mm 
diameter)  at  a  nominal  wire  spacing  of  1.0,  2.0  or  3.0  mm. 
Photographs  of  these  wire-wrapped  tubes  are  shown  in  Figures 
3.9,  3.10,  and  3.11. 
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Figure  3.5    Tracing  of  the  Fin  Profile  of  Tube  54 
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(a)  Rectangular  (Tube  17) 


(b)  Parabolic  (Tube  38) 


(c)  Triangular  (Tube  52) 


(d)  Trapezoidal  (Tube  53) 


Figure  3.7    Cross-Sectional  Photographs  of  Tubes  with 
Different  Fin  Shapes  t^  =  0.5  mm). 
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IV.  SYSTEM  OPERATION  AND  DATA  REDUCTION 

A.   SYSTEM  OPERATION 

Since  copper  has  poor  wetting  characteristics  with 
water,  steam  will  normally  condense  on  copper  under  a 
partial  dropwise  condensation  mode,  which  is  more  effective 
than  the  filmwise  condensation  mode.  Since  the  purpose  of 
this  investigation  was  to  take  data  with  filmwise  condensa- 
tion and  most  of  the  tested  tubes  were  made  from  copper, 
great  care  had  to  be  taken  to  ensure  that  the  filmwise 
condensation  mode  was  in  fact  occurring.  In  order  to  ensure 
this,  the  tubes  had  to  be  treated  according  to  the  following 
procedure : 

1.  The  tube  was  rinsed  with  tap  water  to  remove  any 
contaminants  that  are  soluble  in  the  water. 

2.  A  mixture  of  equal  parts  of  sodium-hydroxide  and  ethyl 
alcohol  was  prepared  and  heated  to  about  80°C,  while 
frequently  being  stirred  until  it  became  watery. 

3.  A  coating  of  this  mixture  was  applied  uniformly 
around  the  tube. 

4.  The  tube  was  placed  in  a  steam  bath  and  was  heated  by 
the  steam  for  about  an  hour. 

5.  A  new  coating  was  applied  to  the  tube  every  10 
minutes . 

6.  The  tube  was  then  rinsed  with  distilled  water  and  put 
immediately  into  the  test  section  to  avoid  any  contam- 
inants depositing  on  the  tube  which  may  lead  to  the 
dropwise  problem. 

This  process  resulted  in  the  formation  of  a  thin  layer  of 
dark  oxide  that  has  high  wetting  characteristics.  Since 
this  layer  was  thin,  its  thermal  resistance  was  negligible. 
This  procedure  was  followed  each  time  prior  to  the  installa- 
tion of  a  tube.  However,  when  re- testing  an  already- 
darkened  tube,  is  was  heated  in  the  steam  only  for  15  to  20 
minutes . 

Following  the  procedure  described  by  Georgiadis  [5]  ,  the 
test   apparatus   was   brought   to   operating   pressure   and 
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temperature  by  adjusting  the  input  power  to  the  boiler 
heaters,  the  cooling  water  flow  rate  through  the  tube,  and 
the  cooling  water  flow  rate  to  the  auxiliary  condenser. 
Steady-state  conditions  were  assumed  when  the  operating 
conditions  were  stabilized  with  a  steam  temperature 
variation  of  ±  2  u.  V  and  a  temperature  rise  of  the  cooling 
water  of  ±  0.005  K  and  ±  0.01  K  for  atmospheric  and  vacuum, 
respectively.  Once  steady-state  operating  conditions  were 
reached,  the  cooling  water  flow  rate  through  the  test  tube 
was  fed  to  the  computer  manually  while  the  temperature  rise 
of  the  cooling  water  through  the  test  tube,  vapor  pressure 
and  temperature  were  gathered  automatically  by  the  data 
aquisition  system.  For  cooling  water  flow  rates  of  80 
percent  (4.44  m/s  for  19  mm  O.D.  and  12.7  mm  I.D.  tubes,  and 
2.84  m/s  for  17.5  mm  O.D.  with  15.6  mm  I.D.  tubes),  70,  62, 
54,  45,-35,  26,  and  20  percent,  and  again  80  percent,  two 
sets  of  data  were  taken.  These  cooling  water  flow  rates 
were  selected  to  give  approximately  equally-spaced  heat  flux 
values.  After  each  change  of  the  cooling  water  flow  rate 
through  the  tube,  the  system  pressure  experienced  a  slow 
drift;  so  an  adjustment  of  the  water  flow  rate  through  the 
auxiliary  condenser  was  required  to  maintain  the  system 
pressure  at  the  operating  pressure.  As  mentioned  in  Chapter 
III,  a  view  port  was  provided  for  visual  observation  to 
ensure  filmwise  condensation.  Before  each  data  collection, 
the  appearance  of  the  film  was  checked.  If  the  film 
appeared  to  be  patchy  or  there  was  an  indication  of  dropwise 
condensation,  the  run  was  discontinued  and  the  data  were 
discarded.  However,  there  were  cases  where  the  film 
appeared  filmwise  but  the  data  collected  at  the  end  of  the 
run  (cooling  water  flow  rate  of  80  percent)  were  different 
from  that  collected  at  the  beginning  for  the  same  flow  rate. 
For  example,  the  heat- transfer  coefficient  was  as  much  as  10 
percent  greater  for  the  last  data  point  and  the  cooling 
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water  temperature  rise  was  also  greater  than  that  measured 
at  the  beginning.  As  discussed  by  Georgiadis  [5],  this 
increased  heat-transfer  coefficient  appears  to  be  a  result 
of  the  tube  undergoing  partial  dropwise  condensationwith 
exposure  to  the  steam.  Since  this  trend  was  observed  even 
though  no  droples  were  visible,  it  is  possible  that  the 
dropwise  condensation  was  taking  place  at  a  microscopic 
level,  especially  near  fin  edges  with  a  very  thin  condensate 
film.  This  phenomenon  was  observed  primarily  for  the  runs 
which  followed  th,e  first  treatment  for  darkening,  probably 
due  to  the  contaminants  of  the  machine  shop  and  since  not  a 
good  layer  of  the  dark  oxide  was  obtained.  All  data 
presented  in  this  thesis  displayed  less  than  3%  disagreement 
in  the  steam-side  heat- transfer  coefficient  between  initial 
and  final  data  sets. 

B.   DATA  REDUCTION 

Initially,  the  program  used  for  data  reduction  was  the 
same  as  that  used  by  Flook  [6]  including  property  functions, 
calibration  curves  for  the  cooling  water  flowmeter  and  for 
all  thermocouples  as  well  as  the  temperature  rise  due  to 
frictional  heating  within  the  mixing  chamber.  Since  tubes 
49,  50  and  51  had  a  different  inside  diameter  than  the  tubes 
tested  previously,  and  tubes  57  thru  62  were  manufactured 
from  metals  with  different  thermal  conductivity,  the  program 
was  modified  to  include  options  for  different  tube  diameter, 
thermal  conductivity,  fin  shapes,  and  Sieder-Tate  constant 
for  the  inside  heat- transfer  coefficient. 

The  separation  of  the  individual  thermal  resistances 
(water-side,  wall,  and  vapor-side)  from  the  overall  heat- 
transfer  resistance  is  very  important  in  order  to  obtain 
expressions  for  the  vapor-side  heat-transfer  coefficient. 
The  overall  heat-transfer  resistance  is  given  by  equation 
(4.1),  while  the  inside  heat-transfer  coefficient  is  given 
by  a  Sieder-Tate  type  equation  (equation  (4.2)).   The  value 
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of  the  leading  coefficient  C^  must  be  known  in  order  to 

calculate  the  inside  heat-transfer  coefficient  from  equation 

(4.2)  and  consequently  the  outside  heat-tranfer  coefficient. 

In  order  to  determine  the  value  of  C_-   two  methods  were 

J-  » 

considered:   the  "direct"  method  and  the  "modified  Wilson 
plot"  method. 
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Li_  (4.1) 


uo.Ao   hi  Ai  '  ho  Ao    Ao 


Nu=iiE£l.ClRe0.8Prl/3[¥]0-",+  B  (4.2) 


3 

1.   "Direct"  Method 

This  method  is  used  to  find  the  leading  coefficient 
for  the  Sieder-Tate  equation  from  an  instrumented  tube.  As 
described  by  Georgiadis  [5]  ,  a  thick-wall  smooth  tube  was 
manufactured  with  six  thermocouples  inserted  into  channels 
around  the  periphery  of  the  tube.  The  average  wall  tempera- 
ture was  found  by  averaging  the  temperatures  indicated  by 
the  six  thermocouples.  He  showed  that  the  wall  temperature 
distribution  followed  a  cosine  curve  given  by  equation  (4.3) 
with  a  maximum  drop  of  18  K  between  the  top  and  bottom  of 
the  tube. 


at       cose  (4.3) 

AT        a  ~ 
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where  a  is  found  to  be  from  0.135  to  0.202  and  from  0.115  to 
0.179  under  vacuum  and  at  atmospheric  pressure,  respec- 
tively, in  order  to  fit  the  temperature  measurements.  Using 
the  average  reading  of  all  the  thermocouples,  Georgiadis  [5] 
found  values  of  0.0635  and  26.4  for  the  leading  Sieder-Tate 
coefficient  and  the  constant  B  in  equation  (4.2),  respec- 
tively, based  on  two  runs  each  under  vacuum  and  at  atmos- 
pheric pressure.  The  value  of  C  (0.0635)  is  greater  than 
the  well-known  Sieder-Tate  constant  of  0.027  for  the  plain 
tubes,  mainly  owing  to  the  coiled  insert.  The  constant 
B  =  26.4  is  used  for  improved  fitting  of  the  experimental 
data. 

2.   "Modified  Wilson  Plot"  Method 

This  method  is  a  modification  of  the  original  Wilson 
plot  method  as  modified  by  Briggs  and  Young  [33]  to  accept 
data  collected  at  various  flow  rates  and  temperatures.  A 
Sieder-Tate  equation  was  used  for  the  inside  heat- transfer 
coefficient,  while  a  Nusselt  type  equation  was  used  for  the 
outside  heat-transfer  coefficient  as  given  by  equation 
(4.4): 


h„   =  3 


kf1  Pf  (Pf  -  Pv)  hfg\  i/3 


°    \    Mf  d0  q 


)   "  (4.4) 


Both  constants  in  equations  (4.2)  and  (4.4)  had  to  be  deter- 
mined iteratively.  Substituting  equation  (4.2)  (with  B 
=  0.0)  and  equation  (4.4)  in  the  equation  for  the  overall 
heat  transfer  resistance  given  by  equation  (4.1),  results  in 
equation  (4.5)  below: 


1  .,  \r.  Es-JL    +  1  (^.s) 


(*. "  "•) 


Ci  kf  Q 
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where 


r  = 


kf    pf    (pf    -    pv)    hfg    |   1/3 

^f    D0   q 


(4.6) 


=  i_Re0.8   Pri/3    [J* 
Di  \Uw 


0.  1<+ 


1/3  (£]  (4'7) 

Equation  (4.5)  is  a  linear  equation  of  the  form: 


Y  = 


mX  +b 

where 


(4.8) 


Y  = 


*,"-).' 


(4.9) 


x  - 


ftp  r 
kf  n 


(4.10) 


m  =  7T 


1 


(4.11) 


1 

b  =  ^r 


(4.12) 


To   begin   the   iteration,   reasonable   values   were 
assumed  for  C^  and  p.   With  these  values,  the  Y  and  X  values 
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were  calculated  and  a  least-square  technique  was  used  to 
compute  the  slope  and  intercept  values  in  equation  (4.8).  A 
new  set  of  values  for  Ci  and  p  was  then  computed  according 
to  equations  (4.11)  and  (4.12).  This  procedure  was  repeated 
until  the  assumed  and  computed  values  for  Ci  and  R  agreed 
to  within  0.1  percent.  Based  on  the  Nusselt  theory,  R 
should  take  a  value  of  0.655,  which  is  true  for  the  case  of 
zero  vapor  shear.  However,  for  the  experimental  conditions 
in  this  thesis,  the  vapor  velocity  was  from  1  m/s  (atmos- 
pheric runs)  to  2  m/s  (vacuum  runs),  thus  resulting  in  R 
values  as  high  as  0.75.  To  account  for  the  vapor  shear 
properly,  a  correlation  developed  by  Fujii  and  Honda  [34] 
was  also  considered  as  shown  below: 


Nu  i/5 

pe  u.b  =  °'96  F  (4.13) 


Using  a  Nusselt-type  equation  for  the  steam-side  coeffi- 
cient, Georgiadis  [5]  found  the  leading  coefficient  (for 
equation  (4.2))  C^  to  be  0.071,  with  the  B  value  set  equal 
to  zero  for  a  smooth  tube.  This  C^  value  resulted  in  an 
inside  heat  transfer  coefficient  up  to  6  percent  greater 
than  that  based  on  the  direct  method.  Flook  later  used  a 
Fujii-type  [34]  equation  instead  of  the  Nusselt-type  equa- 
tion for  the  steam-side  coefficient.  This  resulted  in  a 
slightly  higher  value  (up  to  3  percent)  for  the  leading 
coefficient  Cj.  The  program  used  for  the  data  collection 
allows  an  option  for  selecting  either  the  Fujii-type  or 
Nusselt-type  equation  for  the  steam-side  coefficient. 

Georgiadis  [5]  and  Flook  [6]  thought  that  the  the 
"direct"  method  is  more  reliable,  so  the  values  of  0.0635 
and   26.4   were   used   for   the    constants   C^   and   B, 
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respectively,   for  the  data  reduction  for  tubes  with  an 
inside  diameter  of  12.7  mm. 

Later,  Flook  [6]  used  the  "modified  Wilson  plot" 
method,  with  the  Fujii-type  equation  for  the  steam-side 
coefficient  to  find  the  leading  coefficient  for  the 
Sieder-tate  equation  for  a  copper  tube  and  a  stainless  steel 
tube  with  thin  (i.e.,  0.5  mm  thickness)  tube  walls  since  it 
was  not  possible  to  manufacture  an  instrumented  tube  for 
using  the  "direct"  method  to  find  the  leading  coefficient 
for  the  equation  (4.1).  For  the  copper  and  stainless  steel 
tubes  tube,  values  of  C±  =  0.0756  and  C±  =  0.0688  were 
obtained  respectively.  For  this  thesis  effort,  initially 
the  values  of  C^  =  0.0635  and  B  =  24.6  were  found  with  the 
"direct  method"  as  said  above,  but  finally  the  "Modified 
Wilson  Plot"  directly  on  the  finned  tubes  was  used.  This 
was  the  same  method  as  described  above  but  for  each  tube  a 
different  Sieder-Tate  coefficient  was  found  and  used  for  the 
data  reduction. 
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V.  RESULTS  AND  DISCUSSION 

A.   INTRODUCTION 

During  this  thesis  effort,  a  number  of  data  runs  were 
made  using  the  procedures  described  in  Chapter  IV.  Each 
tube  was  tested  at  least  three  times,  both  under  vacuum  and 
at  atmospheric  pressure  on  different  days,  to  ensure  repeat- 
ability of  the  data.  Complete  filmwise  condensation  condi- 
tions were  maintained,  and  the  non- condensing  gas 
concentration  was  calculated  at  the  beginning  and  at  the  end 
of  each  run  for  every  tube  to  ensure  there  were  no  major 
leaks  in  the  system.  The  computed  mass  concentration  of 
non- condensing  gases  was  kept  between  0  and  -  2.5  percent  in 
order  for  the  data  run  to  be  accepted.  The  mass  concentra- 
tion of  the  non- condensing  gases  was  .  always  a  negative 
number  because  of  the  slight  inaccuracies  in  the  measurement 
of  pressure  and  steam  temperature.  As  discussed  earlier  in 
Chapter  III,  the  test  apparatus  would  allow  only  a  negli- 
gible amount  of  non- condensing  gases  to  be  leaked  into  the 
apparatus.  Since  continuous  venting  was  provided  throughout 
all  runs  (see  Chapter  III),  build  up  of  non- condensing  gases 
was  not  possible. 

A  summary  of  finned  tubes  tested  by  Georgiadis  [5]  and 
Flook  [6]  and  those  tested  during  this  thesis  effort,  as 
well  as  the  resulting  enhancements  are  provided  in  Table 
III.  Further,  Table  IV  presents  the  wire-wrapped  tubes 
tested  and  their  heat- transfer  performance. 
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B.   WATER- SIDE  HEAT-TRANSFER  COEFFICIENTS 

As  mentioned  in  Chapter  IV,  Georgiadis  [5]  used  two 
methods  to  find  the  Sieder-Tate  coefficient:  the  "direct" 
method  and  the  "modified  Wilson  plot"  method  The  "direct" 
method  involved  the  measurement  of  the  average  tube  wall 
temperature  using  six  thermocouples  embedded  within  the  wall 
of  a  smooth  tube.  He  took  data  with  filmwise  condensation 
occurring  outside.  He  changed  the  water  velocity  from  about 
0.8  to  4.5  m/s  and  correlated  the  data  resulting  in  a 
Sieder-Tate  constant  C^  of  0.064  (see  equation  (4.2))  with  a 
B  value  of  26.4.  Also,  taking  data  on  an  unistrumented 
smooth  tube,  and  making  a  modified  Wilson  analysis,  he  found 
a  Sieder-Tate  constant  of  0.071  with  the  B  value  set  equal 
to  zero.  When  the  h^  values  computed  using  the  results  of 
these  two  methods  were  compared,  they  agreed  to  within  6 
percent . 

During  the  present  study,  however,  a  third  approach  was 
tested.  For  this  purpose,  the  modified  Wilson  analysis  was 
carried  out  directly  on  finned-tube  data,  resulting  in  C^ 
values  around  0.069  with  the  B  value  set  equal  to  zero.  The 
h^  values  computed  using  this  analysis  in  fact  lie  between 
the  values  computed  by  the  two  methods  decribed  earlier. 
For  this  reason,  the  third  method  was  used  throughout  this 
investigation  in  computing  the  outside  heat- transfer 
coefficient . 

Since  the  water-side  geometry  for  tubes  45  through  47, 
49  through  51,  52  through  56  ,  57  through  62  and  63  through 
71  was  essentially  the  same,  the  values  of  C^  should  be  very 
nearly  the  same  for  all  these  tubes.  However,  circumferen- 
tial temperature  variations  can  influence  the  inside  coeffi- 
cient, and  these  variations  will  depend  on  the  condensate 
retention  angle.  As  discussed  in  Chapter  II,  the  retention 
angle  is  strictly  dependent  on  the  fin  spacing.  Therefore, 
the  experimentally  found  C-  values  should  vary  from  tube  to 
tube . 
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C.   REPEATABILITY  OF  DATA 

In  order  to  insure  the  reliability  of  the  data  taken, 
all  data  runs  were  repeated,  as  Georgiadis  [5]  and  Flook 
[6]  did,  at  least  three  times  on  different  days.  The 
computed  steam-side  coefficients  for  similar  conditions 
(i.e.,  same  tube  and  about  the  same  operating  conditions)  on 
different  days  agreed  to  within  ±  5  percent  for  some  tubes 
and  ±  10  percent  for  others.  Additionally,  data  runs  were 
performed  on  six  finned  tubes  (tubes  6,  17,  27,  28,  36,  38) 
under  similar  conditions  to  verify  the  repeatability  with 
data  taken  by  Flook  [6]  .  Georgiadis  [5]  also  had  tested 
tubes  6  and  17.  Georgiadis  and  Flook  processed  their  data 
using  the  Sieder-Tate  constant  found  by  the  "direct"  method 
(see  Section  B  of  Chapter  IV).  Therefore,  in  order  to 
perform. a  fair  comparison,  their  data  were  reprocessed  by 
the  method  used  during  this  investigation  (i.e.,  using  the 
"modified  Wilson  plot"  method  directly  on  finned  tube  data). 
Figure  5.1  shows  the  experimental  steam-side  heat-transfer 
coefficients  of  Georgiadis  [5]  and  Flook  [6]  and  those 
obtained  during  this  investigation  for  tube  number  6  under 
vacuum.  Figure  5.2  shows  similar  data  for  tube  6  under 
atmospheric  pressure,  whereas  Figures  5.3  and  5.4  show 
similar  data  for  tube  17.  For  comparison  purposes,  the 
smooth  tube  data  were  plotted  as  well  as  the  data  predicted 
by  Nusselt  theory.  As  can  be  seen  from  Figure  5.1,  the  data 
obtained  during  this  work  under  vacuum  conditions  fall  above 
and  below  the  data  of  Georgiadis  [5]  and  Flook  [6]  with  a 
maximum  variation  of  20  percent.  At  atmospheric  conditions, 
the  data  fall  about  15-20  percent  below  those  of  Georgiadis 
[5]  .  The  agreement  with  tube  17  is  much  better  as  can  be 
seen  from  Figures  5.3  and  5.4.  For  tubes  27,  28,  36,  and  38 
(see  Figures  5.5  through  5.9)  the  agreement  is  very  good, 
with  a  deviation  of  only  ±  5  percent  except  for  tube  number 
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36  at  atmospheric  pressure.  In  this  case,  the  agreement  was 
±  26  percent.  The  disagreement,  especially  for  tube  6,  is 
probably  due  to  a  very  small  contamination  on  the  outside 
tube  surface  area  of  the  tube  leading  to  partial  dropwise 
condensation  conditions  during  the  runs  made  by  Georgiadis 
and  Flook.  Typical  uncertainty  bands  are  also  included  at 
low  and  high  heat  flux  to  indicate  the  maximum  possible 
uncertainty.  However,  the  repeatability  for  most  of  the 
tubes  shows  that  the  experimental  uncertainty  is  always  less 
than  that  indicated  by  the  uncertainty  bands. 

In  these  figures,  the  steam-side  coefficient  is  plotted 
versus  the  heat  flux,  and  as  is  always  the  case  for  conden- 
sation, the  heat-transfer  coefficient  decreases  as  the  heat 
flux  increases.  The  curves  shown  in  these  figures  (and 
subsequent  figures)  are  the  least-squares-fit  curves 
according  to  the  following  equation: 


q  =  a  AT^  (5.1) 

■ 

where  AT  was  computed  using  the  following  equation: 


q  =  h0  AT  (5.2) 

where  q  was  measured  experimentally  and  the  steam-side  coef- 
ficient was  calculated  by  subtracting  the  inside  and  wall 
thermal  resistances  from  the  overall  resistance  as  given  by 
equation  (5.3): 


where 


h°  =   I   ™    5Ti  (5-3) 

U"  *«  "  Ai  hi 
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Rw  - 


DQ  lnCDp/Dj) 
2  kw 


(5.4) 


In  equation  (5.1)  the  coefficient  a  and  the  exponent  b  are 
experimentally  determined  constants.  The  values  of  a  and  b 
both  under  vacuum  and  atmospheric  conditions  for  all  the 
finned  tubes  and  for  the  wire-wrapped  tubes  tested  are  given 
in  Tables  V  and  VI,  respectively. 

TABLE  V 
CONSTANTS  OF  EQUATION  (5.1)  FOR  FINNED  TUBES  TESTED 


Tube 

VACUUM   RUNS 

ATMOSPHERIC  RUNS 

Tube 
No. 

s 
(mm) 

c 
(mm) 

e 
^  mm ) 

a 

b 

a 

b 

06 

1.50 

1.00 

1.00 

61190 

0.72588 

91685 

0.71222 

17 

1.50 

0.50 

1.00 

63210 

0.72465 

95079 

0.71433 

27 

1.50 

1.00 

0.50 

49315 

0.71454 



28 

2.00 

1.00 

0.50 

52515 

0.71739 



36 

0.00 

2.10 

1.00 

63847 

0.72346 

93246 

0.71243 

38 

1.50 

0.50 

1.00 

67381 

0.72661 

102870 

0.71715 

45 

0.00 

2.50 

1.00 

63459 

0.72502 

99843 

0.70538 

46 

0.00 

1.60 

1.00 

65976 

0.72740 

101200 

0.71728 

47 

0.00 

1.06 

1.00 

63456 

0.72683 

98200 

0.71167 

49 

0.51 

0.34 

1.00 

61325 

0.72920 

89327 

0.71046 

50 

28144 

0.70192 

43092 

0.64361 

51 

0.59 

0.32 

1.00 

61928 

0.72950 

93058 

0.71267 

52 

1.50 

0.50 

1.00 

63672 

0.72548 

95028 

0.71388 

53 

1.50 

0.50 

1.00 

62772 

0.72613 

94397 

0.71270 

54 

1.50 

1.00 

1.00 

61452 

0.72125 

94299 

0.71061 

55 

1.50 

1.00 

1.00 

66987 

0.72326 

90906 

0.  70660 

56 

1.50 

1.00 

1  .00 

65674 

0.72303 

94829 

0.70866 

57 

0.00 

2.  10 

1.00 

65871 

0.72540 

94095 

0.69296 

58 

0.00 

2.  10 

1.00 

42916 

0.70566 

73831 

0.68875 

59 

0.00 

2.  10 

1.00 

24364 

0.67314 

38902 

0.67594 

60 

1.50 

1.00 

1.00 

57695 

0.71805 

82971 

0.69082 

61 

1.50 

1.00 

1.00 

40173 

0.70228 

72273 

0.69034 

62 

1.50 

1.00 

1.00 

50021 

0.71093 

85279 

0.69645 

D.   EFFECT  OF  FIN  PITCH  ON  HEAT  TRANSFER  PERFORMANCE  OF 
SPIRAL  TUBES  WITH  TRIANGULAR- SHAPED  FINS 

This  section  presents  results  showing  the  variations  of 

the  steam-side  heat- transfer  coefficient  with  heat   flux 

having  fin  pitch  as  a  parameter.   Data  were  taken  on  four 
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TABLE  VI 

CONSTANTS  OF  EQUATION  (5.1) 
FOR  WIRE-WRAPPED  TUBES  TESTED 


Tube 

VACUUM  RUNS 

ATMOSPHERIC  RUNS 

Tube 
No. 

Dw 
(mm) 

s 
(mm) 

(mm) 

a 

b 

a 

b 

63 
64 
65 
66 
67 
68 
69 
70 
71 

1.6 
1.6 
1.6 
1.0 
1.0 
1.0 
0.5 
0.5 
0.5 

0.94 
2.03 
3.02 
0.95 
1.82 
2.91 
1.10 
2.04 
3.13 

2.54 
3.63 
4.62 
1.95 
2.82 
3.91 
1.60 
2.54 
3.63 

26207 
32515 
32815 
27097 
35403 
36521 
34291 
40254 
39179 

0.68301 
0.69452 
0.69733 
0.68988 
0.70265 
0.70144 
0.69228 
0.70477 
0.70190 

46280 
52572 
52403 
47762 
58145 
58421 
59000 
63569 
62272 

0.65983 
0.67158 
0.66672 
0.65921 
0.67488 
0.67219 
0.66327 
0.67866 
0.67547 

copper  tubes  with  spiral  triangular  fins.  These  tubes  have 
the  same  fin  height  of  1.0  mm  and  have  fin  pitches  of  1.06, 
1.6,  2.1,  and  2.5  mm.  Figures  5.10  and  5.11  present  data 
for  these  tubes  under  vacuum  and  at  atmospheric  pressure, 
respectively.  The  smooth-tube  data  and  a  curve  representing 
Nusselt  theory  also  included  for  comparison.  The  best  heat- 
transfer  performance  was  obtained  with  the  tube  with  a  fin 
pitch  of  1.6  mm. 

As  shown  in  Appendix  C,  the  uncertainty  in  the  calcula- 
tion of  the  steam-side  coefficient  increases  as  the  heat 
flux  decreases.  Therefore,  the  comparison  of  the  finned 
tubes  should  be  performed  at  a  high  heat  flux,  where  the 
uncertainty  is  small.  The  comparison  of  finned  tubes  is 
made  through  the  enhancement  ratio,  Eo .  This  ratio  is 
defined  as  the  steam-side  heat- transfer  coefficient  of  a 
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finned  tube  to  that  of  the  smooth  tube  (same  diameter  as  the 
finned  tube  root  diameter)  at  the  same  heat  flux.  Heat  flux 
values  of  0.25  and  0.75  MW/m  were  chosen  for  vacuum  and  at 
atmospheric  conditions,  respectively.  For  the  spiral  trian- 
gular fins,  maximum  enhancement  ratios  of  about  3.9  and  6.1 
under  vacuum  and  at  atmospheric  pressure  were  found.  The 
enhancement  ratio  at  atmosheric  pressure  is  always  higher 
than  that  under  vacuum  conditions.  At  atmospheric  pressure, 
a  higher  temperature  exists,  so  the  condensate  has  a  smaller 
viscosity,  which  results  in  improved  drainage  from  the  fin 
valleys,  and  smaller  surface  tension  which  results  in  a 
smaller  retention  angle.  As  discussed  in  Chapter  II,  the 
flooded  portion  of  the  tube  has  another  thermal  resistance 
due  to  the  thick  layer  of  condensate.  Reducing  the  flooded 
portion  of  the  tube  increases  the  heat-transfer  performance. 
Cross  plots  of  the  enhancement  ratio  versus  fin  pitch 
are  shown  in  Figure  5.12,  while  Figure  5 .  13  shows  a  cross 
plot  of  the  normalized  ratio  Eo/Ar  (the  ratio  of  the 
enhancement  ratio  to  the  area  ratio).  Generally,  as  the  fin 
pitch  increases,  Eo/Ar  increases.  Table  III  shows  that  tube 
47,  with  a  fin  pitch  of  1.06  mm,  has  the  largest  area  ratio, 
while  tube  45  with  fin  pitch  of  2.5  mm  has  the  smallest  area 
ratio,  and  tube  47  has  a  poorer  performance  than  either 
tubes  45  or  46.  The  poor  performance  shown  by  tube  47  can 
be  explained  by  the  fact  that,  as  the  pitch  decreases,  the 
area  of  the  tube  increases,  but  at  the  same  time  the  reten- 
tion angle  increases  and  more  flooding  occurs.  This  means 
that  as  the  fin  pitch  increases,  the  interfin  spacing  was 
covered  by  a  thick  layer  of  condensate.  The  additional 
thermal  resistance  induced  by  this  layer  of  condensate  over- 
powers the  benefit  gained  from  increased  surface  area,  so 
the  heat-transfer  performance  is  reduced. 

As  the  pitch  increases  from  1.05  mm  to  1.6  mm,   the 
retention  angle  decreases  more  than  the  area  decreases. 
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This  results  in  a  larger  enhancement  ratio  for  the  tube  with 
a  pitch  of  1.6  mm.  Beyond  this  point,  the  area  ratio 
decreases  while  the  retention  angle  decreases  slowly  and 
results  in  a  smaller  heat  transfer  performance.  In  order  to 
obtain  a  clearer  understanding  on  the  heat-transfer  perform- 
ance, the  enhancement  ratio  was  divided  by  the  area  ratio, 
and  thus,  the  effect  of  the  changing  area  was  eliminated. 
As  Figure  5.13  shows,  the  enhancement  ratio  depends  on  other 
factors  in  addition  to  the  fin  area,  such  as  the  surface 
tension  effect.  As  discussed  in  Chapter  II,  the  surface 
tension- induced  pressure  gradient  from  the  fin  tip  to  the 
fin  root  is  responsible  for  thinning  of  the  condensate  film 
and  thereby  improving  the  heat-transfer  performance  in  the 
unflooded  portion  of  the  tube.  Also,  the  surface-tension 
forces  cause  condensate  flooding,  resulting  in  poorer 
performance  in  the  flooded  portion  of  the  tube.  Figure  5.13 
shows  that  the  normalized  enhancement  ratio  is  higher  for 
the  tube  with  fin  pitches  of  2.1  mm  and  2.5  mm  under  vacuum 
and  at  atmospheric  pressure,  respectively.  Therefore,  the 
optimum  fin  pitch  is  between  2.1  and  2.5  mm  based  on  normal- 
ized enhancement  ratio,  while  the  optimum  fin  pitch  is 
1.6  mm  based  on  the  enhancement  ratio.  As  shown  by  Edwards 
et  al.  [20],  as  the  pitch  increases  for  the  same  fin  height, 
the  heat- transfer  coefficient  increases.  However,  as  the 
fin  pitch  increases,  the  tube  is  easily  flooded.  Because  the 
retention  angle  is  greater  under  vacuum  than  that  at  atmos- 
pheric pressure,  the  tube  with  fin  pitch  of  2 . 5  mm  has  a 
smaller  retention  angle  than  that  of  the  tube  with  fin  pitch 
of  2.1  mm  under  vacuum,  resulting  in  poorer  heat- transfer 
performance.  However,  at  atmospheric  pressure,  since  the 
retention  angle  is  less  than  under  vacuum  conditions,  the 
tube  with  a  fin  pitch  of  2.5  mm  has  a  better  heat- transfer 
performance  than  the  tube  with  a  fin  pitch  of  2.1  mm. 
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E.   EFFECT  OF  FIN  SHAPE  ON  HEAT-TRANSFER  PERFORMANCE 

In  order  to  study  the  effect  of  fin  shape,  data  were 
taken  on  two  sets  of  copper  tubes  with  fins  of  four 
different  shapes.  All  the  fins  were  manufactured  with  same 
fin  spacing  and  fin  height.  The  first  set  of  tubes  (17,  38, 
52  and  53)  had  rectangular,  parabolic,  triangular  and  trape- 
zoidal fin  shapes,  respectively,  with  a  fin-base  thickness 
of  0.5  mm,  while  the  second  set  of  tubes  (6,  54,  55  and  56) 
had  a  fin-base  thickness  of  1.0  mm.  Dimensions  for  these 
fins  are  given  in  Table  III. 

The  performance  of  tubes  17,  38,  52,  and  53  under  vacuum 
conditions  is  shown  in  Figure  5.14,  while  Figure  5.15 
depicts  their  performance  at  atmospheric  pressure.  For 
comparison  purposes,  data  for  a  smooth  tube  are  also  shown. 
The  tube  with  the  "parabolic"  fin  profile  (tube  38)  showed 
the  best  heat- transfer  performance,  while  the  other  three 
tubes  performed  about  equally,  under  both  pressure  condi- 
tions. As  shown  in  Table  III,  an  enhancement  ratio  of  4.1 
and  6.2  were  obtained  for  the  tube  with  "parabolic"  fins 
under  vacuum  and  at  atmospheric  pressure,  respectively. 
Also,  it  can  be  seen  that  the  area  ratio  of  the  tube  with 
"parabolic"  fins  is  less  than  that  of  the  tube  with  rectan- 
gular fins  and  the  tube  with  trapezoidal  fins,  and  larger 
than  that  of  the  tube  with  triangular  fins.  Eliminating  the 
effect  of  increased  area,  the  ratio  Eo/Ar  is  larger  for  the 
tube  with  "parabolic"  fins  under  atmospheric  and  vacuum 
conditions . 

The  reason  for  the  greater  enhancenent  ratio  is  probably 
due  to  the  continuous  change  of  radius  of  curvature 
(increasing  from  the  fin  tip  to  the  fin  root)  for  the  "para- 
bolic" shaped  fins.  The  condensate  film  has  a  convex  shape 
at  the  fin  tip  and  a  concave  shape  at  the  fin  root.  The 
condensate  film  follows  approximately  the  curvature  of  the 
wall  surface  at  the  fin  tip.  Because  of  the  convex  conden- 
sate surface  at  the  fin  tip  and  the  concave  condensate 
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surface  near  the  fin  root,  the  pressure  within  the  conden- 
sate is  larger  at  the  fin  tip  and  smaller  at  the  fin  root. 
Therefore,  an  appreciable  pressure  gradient  exists  from  the 
fin  tip  to  the  fin  root.  The  gradual  increase  of  radius  of 
curvature  results  in  a  gradual  decrease  in  pressure  within 
the  condensate,  which  is  very  important  for  improved  conden- 
sate flow,  resulting  in  a  thinner  film  and  larger  heat- 
transfer  coefficient  than  if  the  fin  sides  were  flat. 
Therefore,  the  parabolic  fins  should  outperform  all  other 
three  tubes,  as  shown  in  Figures  5.14  and  5.15. 

Also,  Adamek  [24]  and  Mori  et  al.  [8]  have  shown  that 
the  optimum  fin  shape  is  that  which  induces  a  continuous 
pres'sure  gradient  due  to  the  surface  tension  effect  or  which 
has  large  curvature  at  the  fin  tip  and  continuously- 
decreasing  toward  the  fin  root.  Therefore  a  continous 
decrease  in  the  pressure  gradient  exists  and  this  thins  the 
condensate  film  continuously,  resulting  in  better  heat 
transfer  performance.  For  the  case  of  tube  38,  there  is  a 
continuous  decrease  of  the  curvature,  while  this  was  not 
happening  for  the  other  three  tubes  of  the  first  set 
although  they  had  a  sharp  leading  edge.  Therefore  tube  38 
exhibits  better  heat-transfer  performance  than  the  other 
three  tubes  in  its  group. 

Figures  5.16  and  5.17  show  the  performance  of  tubes  6, 
54,  55  and  56  under  vacuum  and  at  atmospheric  conditions, 
respectively.  Again,  the  smooth  tube  data  are  included  for 
comparison  purposes.  For  this  set  of  tubes,  the  best 
performance  is  obtained  from  tube  55  with  the  trapezoidal 
fins,  while  the  tube  with  triangular  fins  performed  second 
and  the  remaining  two  tubes  performed  about  equally  well 
under  vacuum  conditions.  However,  the  parabolic  and  trian- 
gular shapes  outperformed  the  other  shapes  for  atmospheric 
conditions.  The  poorer  performance  of  tube  54  under  vacuum 
conditions  was  not  expected.   As  Figure  3.7  shows,  the  fins 
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do  not  have  the  shapes  as  claimed  above  due  to  the  difficul- 
ties encountered  in  machining.  The  unexplainable  trends 
shown  by  the  tube  with  "parabolic"  fins  (tube  54)  was  found 
to  be  the  result  of  the  actual  fin  shape  that  was  very 
different  from  what  was  expected.  At  the  conlusions  of  the 
data  runs  presented  in  this  thesis,  this  tube  was  destroyed 
and  a  magnified  photogragh  of  the  fin  cross- section  was 
taken.  As  can  be  seen  from  Figure  3.5,  these  fins  do  not 
have  a  profile  with  gradually  decreasing  curvature  from  the 
fin  tip  to  root.  Therefore,  unlike  in  the  previous  set  of 
tubes  (17,  38,  52  and  53),  the  data  taken  on  the  second  set 
of  tubes  are  inconclusive.  Since  the  condensation  process 
on  a  finned  tube  is  extremely  complex,  owing  to  the  very 
large  number  of  parameters,  it  may  be  unwise  to  draw  conclu- 
sions from  the  above-mentioned  results.  Nevertheless,  the 
data  for  the  second  set  of  tubes  are  presented  in  this 
thesis  for  completeness. 

From  Figure  3.7  it  is  clear  that:  1)  the  "parabolic" 
fins  of  the  tube  38  with  fin  base  thickness  of  0.5  mm  had  a 
straight  fin  side  and  a  fin  tip  with  a  sharp  leading  edge, 
while  careful  examination  of  a  cross-section  of  the  "para- 
bolic" fins  of  tube  54  shows  a  straight  fin  side  with  a 
near  semicircle  at  the  fin  tip.  Tube  38  had  a  small  radius 
of  curvature  at  the  fin  tip,  while  tube  54  had  larger 
radius  of  curvature  at  the  fin  tip.  These  differences  in 
geometry  may  have  caused  the  observed  data. 

F.   EFFECT  OF  FIN  THERMAL  CONDUCTIVITY  ON  PERFORMANCE 

To  investigate  the  effect  of  fin-metal  thermal  conduc- 
tivity on  the  heat- transfer  performance,  four  spirally 
finned  tubes  with  triangular  fins  and  two  tubes  with  rectan- 
gular fins  were  manufactured.  As  shown  in  Table  I,  the 
four  spiral  tubes  were  made  of  copper,  copper-nickel, 
stainless  steel,  and  aluminum,  respectively  (tubes  57,  58, 
59  and  60).  The  two  tubes  with  rectangular  fins  were  made 
of  copper-nickel  and  aluminum  (tubes  61  and  62).   All  tubes 
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had  the  same  fin  height  of  1.0  mm.  The  tubes  with  a  rectan- 
gular fin  shape  have  the  same  fin  dimensions  as  the 
"optimum"  copper  tube  (tube  6)  found  by  Georgiadis  [5], 
while  the  spiral  tubes  had  a  fin  pitch  of  2.1  mm.  The 
results  for  data  runs  taken  under  vacuum  and  at  atmospheric 
pressure  are  shown  in  Figures  5.18  and  5.19,  respectively, 
for  the  spiral  tubes.  For  the  tubes  with  rectangular  fins, 
the  variation  of  heat-transfer  coefficient  with  heat  flux  is 
shown  in  Figures  5.20  and  5.21  under  vacuum  and  at  atmos- 
pheric conditions,  respectively.  Figure  5.18  shows  that  the 
copper  tube  exhibits  the  best  heat- transfer  performance, 
while  the  stainless  steel  shows  the  worst  performance.  The 
second  best  is  the  aluminum  spiral  tube  followed  by  the 
copper-nickel  tube.  The  same  trend  is  also  seen  in  Figure 
5.19  at  atmospheric  pressure.  Enhancement  ratios  as  high 
as  3.5  and  4.4  under  vacuum  and  atmospheric  conditions, 
respectively,  were  found.  Cross  plots  of  enhancement  ratio 
Eo  and  normalized  enhancement  ratio  Eo/Ar  versus  the  thermal 
conductivity  are  shown  in  Figures  5.22  and  5.23, 
respectively. 

A  similar  trend  exists  for  tubes  with  rectangular  fin 
profiles  as  shown  in  Figures  5.20  and  5.21.  For  comparison, 
the  data  of  Flook  [6]  for  the  copper  tube  with  DQ  =  14.5  mm 
and  D^  =  13.5  mm  with  the  same  fin  dimensions  as  the 
aluminum  and  copper-nickel  tube  (tube  39)  are  also  included. 
Also  in  the  same  Figures  the  data  of  tube  6  are  included. 

The  thermal  conductivity  of  aluminum  (167  W/m.K)  is 
about  half  of  copper  (385  W/m.K),  while  copper  nickel  and 
stainless  steel  have  much  lower  values  (i.e.,  45  W/m.K  and 
15  W/m.K,  respectively).  Since  the  thermal  resistance 
through  the  fin  increases  (i.e.,  the  fin  efficiency 
decreases)  with  decreasing  thermal  conductivity,  the  copper 
tube  must  show  the  best  heat- transfer  performance,  while  the 
stainless  steel  tube  must  show  the  poorest  performance.   As 
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can  be  seen  from  Figures  5.18  through  5.22,  this  trend  is 
very  clear. 

G.   PERFORMANCE  OF  WIRE  WRAPPED  TUBES 

1.   Condensate  Retention  Angle  for  Wire-Wrapped  Tubes 

Wire-wrapped  tubes  are  somewhat  similar  to  finned 
tubes  with  regard  to  their  susceptibility  to  condensate 
retention.  Since  the  portion  of  the  tube  with  the  retained 
condensate  would  result  in  a  poor  heat- transfer  performance, 
it  was  necessary  to  study  this  phenomenon  on  these  wire- 
wrapped  tubes.  For  this  purpose,  it  is  possible  to  use  the 
Webb  et  al.  [17]  model  (as  discussed  in  Chapter  II),  which 
was  developed  for  a  tube  with  fins  of  arbitrary  shape  (see 
equation  (2.10)).  However,  they  did  not  test  this  model  for 
wire-wrapped  tubes.  Therefore,  it  was  necessary  to  experi- 
mentally measure  the  condensate  retention  angle  for  these 
tubes,  so  that  the  Webb  et  al.  model  can  be  modified  to 
predict  it  for  wire-wrapped  tubes.  For  this  purpose,  meas- 
urement of  the  condensate  retention  angles  were  made  for  the 
wire-wrapped  tubes  under  static  conditions.  Photographs 
were  taken;  slides  were  made  and  accurate  measurements  were 
made  on  the  screen.  Figure  5.24  shows  photographs  of 
portions  of  tubes  with  a  wire  diameter  of  0.5  mm  and  pitches 
of  1.6,  2.5  and  3.6  mm  (tubes  69,  70  and  71),  respectively. 
Careful  examination  of  the  photographs  revealed  that  water 
was  trapped  all  around  the  tube  in  a  region  very  near  of  the 
wires.  The  film  thickness  was  about  the  same  along  the 
distance  from  the  top  of  the  tube  to  the  bottom,  until  some 
point  where  the  meniscus  of  the  water  had  a  concave  profile 
in  the  region  where  the  tube  starts  to  be  fully  flooded 
between  two  adjacent  wires.  In  this  case,  the  retention 
angle  is  defined  as  the  angle  from  the  bottom  of  the  tube  to 
point  A  (see  Figure  5.24),  where  the  surface  tension  forces 
balance  the  weight  of  the  condensate.  The  retention  angle 
was  computed  by  equation  (5.5)  based  on  the  height  from  the 
bottom  of  the  tube  to  point  A  and  the  tube  radius. 
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Tube  # 

# 

63 

64 

65 

66 

67 

68 

69 

70 

71 

TABLE  VII 
MEASURED  RETENTION  ANGLES  (4>) 


Measured 
Degrees 

116 

81 

73 
118 

87 

70 
102 

68 

47 


Dw 

/   Pv 

(mm) 

(mm) 

1.6 

2.5 

1.6 

3.6 

1.6 

4.6 

1.0 

2.0 

1.0 

2.8 

1.0 

3.9 

0.5 

1.6 

0.5 

2.5 

0.5 

3.6 

Ca 

lculated 

De 

grees 

117 

72 

57 

123 

78 

59 

110 

71 

54 

♦  -  (1  -  |>  (5-5) 


Table  VII  lists  the  results  for  all  nine  tubes  (tubes  63 
through  70).  Using  the  Webb  et  al.  model  [14],  attempts 
were  made  to  predict  the  measured  retention  angles,  and 
comparison  is  shown  in  Figure  5.25  (see  the  "triangular" 
symbols--the  "star"  symbols  will  be  discussed  below).  As 
can  be  seen,  the  Webb  et  al.  model  overpredicts  the  conden- 
sate retention  angle  for  most  of  the  tubes.  Therefore,  this 
model  was  modified  for  the  present  study  (i.e.,  steam 
condensation  on  horizontal  wire-wrapped  tubes).  The  exami- 
nation of  equation  (2.10),  reveals  that  the  condensate 
retention  angle  would  decrease  with  increasing  wire  diameter 
while  all  other  parameters  are  kept  constant.  However,  the 
experimentally  measured  retention  angles  shown  in  Table  VII 
show  the  opposite  trend  (i.e.,  \y  increases  with  increasing 
D  ).   Therefore,  it  was  necessary  to  modify  the  Webb  et  al. 


model  as  shown  in  equation  (5.6) 


I",  -  .      2  °(F'-^ (2»\  «"]  (5.6) 

L       no  Pf  8  t(tb  +  s)  e  -  Ap]  \s    j 


ty    ■  COS  l 


where  a  and  (3  are  empirical  constants  to  fit  the  measured 
retention  angle.  In  order  to  compute  these  a  and  b  values, 
a  numerical  procedure  was  followed,  by  minimizing  the  sum  of 
squares  of  the  deviations  of  the  computed  and  measured  y 
values.  This  procedure  resulted  in  a  and  j3  values  of  0.85 
and  0.18,  respectively.  The  comparison  between  the  values 
computed  by  equation  (5.6)  and  the  experimental  values  is 
shown  in  Figure  5.25.  As  can  be  seen,  this  equation  agrees 
to  within  ±  15  percent  with  the  experimental  data. 
2 .   Experimental  Data  of  Wire-Wrapped  Tubes 

This  section  presents  results  showing  the  variation 
of  the  steam-side  coefficient  with  heat  flux  for  the  wire- 
wrapped  tubes  tested.  Three  wire  diameters  (0.5,  1.0  and 
1.6  mm)  were  used  with  three  different  pitches  for  each  wire 
diameter.  These  are  shown  in  Table  IV;  The  tubes  with  a 
1.6  mm  wire  diameter  had  pitches  of  2.5,  3.6,  and  4.6  mm; 
the  tubes  with  1.0  mm  wire  diameter  had  pitches  of  2.0,  2.8, 
and  3.9  mm;  and,  the  tubes  with  0.5  mm  wire  diameter  had 
pitches  of  1.6,  2.5,  and  3.6  mm.  Figures  5.26  and  5.27  show 
the  variation  of  heat-transfer  coefficient  with  heat  flux 
under  vacuum  and  at  atmospheric  pressure,  respectively,  for 
the  tubes  wrapped  with  a  1 . 6-mm-diameter  wire,  while  Figures 
5.28  and  5.29  show  the  variation  of  the  heat- transfer  coef- 
ficient under  vacuum  and  at  atmospheric  conditions,  respec- 
tively, for  the  tubes  wrapped  with  1.0-mm-wire  diameter. 
Figures  5.30  and  5.31,  show  similar  results  for  the  tubes 
wrapped  with  0 . 5-mm-diameter  wire  under  vacuum  and  at  atmos- 
pheric pressure,  respectively.  Figures  5.26  to  5.31  show 
that  the  best  performance  was  obtained  for  the  tube  with  a 
wire  diameter  of  0.5  mm  and  a  pitch  of  2.5  mm,  both  under 
vacuum  and  at  atmospheric  pressure.  The  second  best  perform- 
ance was  given  by  the  tube  with  a  wire  diameter  of  0.5  mm 
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and  a  pitch  of  3.6  mm  under  both  pressure  conditions.  The 
third  best  performance  was  given  by  the  tube  with  a  wire 
diameter  of  1.0  mm  and  a  pitch  of  2.8  mm  under  atmospheric 
pressure,  and  the  tube  with  a  wire  diameter  of  1.0  mm  and  a 
pitch  of  3.9  mm,  under  vacuum  pressure.  The  tubes  with  wire 
diameters  of  1.6  mm  and  1.0  mm  with  pitches  of  2.5  and 
2.0  mm,  respectively,  showed  the  worst  performance  among  all 
the  tubes  tested. 

As  Figures  5.26  and  5.27  show,  the  performance  of 
the  tubes  with  a  wire  diameter  of  1.6  mm  and  pitches  of  3.6 
and  4.6  mm  interchanged  as  the  pressure  conditions  were 
changed  from  vacuum  to  near- atmospheric.  The  same  trend 
happened  for  the  tube  with  wire  diameter  of  1.0  mm  and 
pitches  of  2.8,  and  3.9  mm.  This  behavior  can  only  be 
explained  by  the  retention  angle  phenomenon.  The  tubes  were 
wrapped  in  order  to  take  advantage  of  the  condensate  thin- 
ning as  a  result  of  the  varying  surface-tension  forces  in 
the  space  between  the  wires.  However,  as  mentioned  in 
Chapter  II,  these  surface-tension  forces  lead  to  a  deleter- 
ious effect  owing  to  condensate  retention,  especially  when 
the  fin  spacing  is  small.  If  the  extent  of  condensate 
retention  is  large,  the  enhancement  gained  over  the 
unflooded  portion  of  the  tube  (owing  to  condensate  thinning) 
may  be  offset  in  poor  performance  in  the  flooded  portion 
with  retained  condensate.  As  discussed  in  subsection  1 
above,  the  retention  angle  is  higher  under  vacuum  conditions 
for  tubes  with  a  pitch  of  3.6  than  for  tubes  with  a  pitch  of 
4.6.  This  is  true  for  tubes  with  pitches  of  2.8  and  3.9  mm 
as  well.  However,  under  atmospheric  pressure,  the  retention 
angle  decreases  resulting  in  less  flooding,  and  better 
thermal  performance.  Figures  5.26  and  5.28  show  that  the 
tubes  with  wire  diameters  of  1.6  and  1.0  mm  with  pitches  of 
2.5,  and  2.0  mm,  respectively,  had  worse  performance  than 
the  smooth  tube  under  vacuum.   Again  these  are  the  results 
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of  the  surface  tension  effects.  For  these  tubes  the  reten- 
tion angle  was  about  100  degrees,  resulting  in  about  half 
of  the  tube  being  effective  for  the  heat  transfer.  The 
other  half  contributes  a  small  amount  to  the  heat  transfer 
performance  due  to  the  thick  layer  of  the  condensate.  The 
heat-transfer  enhancement  ratio  is  shown  clearly  in  Figure 
5.32,  as  a  function  of  p/dw.  The  optimum  p/Dw  appears  to  be 
near  5  to  6. 

3 .   Modifications  to  Fujii  et  al.  Model 

As  mentioned  in  Chapter  II,  Fujii  et  al .  [32]  devel- 
oped a  semi- theoretical  expression  to  predict  succesfully 
their  condensation  data  on  wire-wrapped  tubes,  using  ethanol 
and  R-ll  as  the  working  fluids.  As  stated  in  Chapter  II 
(Section  E),  the  vapor-side  enhancement  could  be  easily 
computed  using  equation  (2.65)  provided  A  >  15 .  However, 
for  most  of  the  tubes  tested  during  this  investigation,  the 
values  of  A  were  as  low  as  6.  Thus,  it  was  not  possible  to 
use  equation  (2.62)  and  equation  (2.59)  had  to  be  used. 
When  equation  (2.59)  is  plotted  as  a  function  of  A  and  q>, 
the  result  is  shown  in  Figure  (5.33).  Notice  that  for 
A  >  18 ,  F-^((p,A)  is  independent  of  cp  Also  for  A  >  18  the 
functional  dependence  of  F^  on  the  angle  <p  can  be  approxi- 
mated by  a  straight  line: 


FxU.A)  «  FxCO.A)  -  m  4.  (5.7) 


where  m  is  the  slope  of  the  straight  line.  In  equation 
(5.7),  Fi(0,A)  for  all  values  of  A  is  given  by  equation 
(2.62). 

Because  of  the  flooding  that  occurs  on  the  lower 
portion  of   the  wire-wrapped   tubes  when  steam   is   being 
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condensed,  in  order  to  compute  the  heat-transfer  performance 
of  the  unflooded  portion  of  the  tube,  the  function  F]((p,A) 
must  be  integrated  from  <p  =  0.0  to  <p  =  k  -  \y,  as  needed  for 
equation  (5.8): 

F2(*.A>  "  f  o/*  Fl(t>A)  d*  (5.8) 

Notice  that  this  equation  is  the  same  as  equation  (2.58) 
except  that  the  integration  is  performed  up  to  the  angle  q> 
(which  is  n  -  \j/).  Fujii  et  al.  neglected  the  retention- 
angle  effect  for  their  low-surface-tension  fluids.  As 
mentioned  earlier,  for  A  >  15,  equation(5 . 7 )  is  valid,  but 
the  slope  m  can  depend  on  the  value  of  A.  In  order  to  find 
the  dependence,  the  slope  was  approximated  for  values  of 
A  >  5  and  these  slopes  are  plotted  in  Figure  (5.34)  as  a 
function  of  A.  Using  a  least  square  fit  of  the  calculated 
data,  the  following  functional  form  was  derived: 


b        c  (5.9) 

m  =  a  A     ,  A  >  5 


where    a  =   0.177,  and 

b  =  -0.756 
Notice  that  the  actual  computation  of  m  was  performed  only 
for  5  <  A  <  19.  The  numerical  integration  of  F^(<p,A)  for  A 
>  18  was  not  possible  owing  to  overflow  limitations  of  the 
computer.  Further,  the  integration  of  Fi(q),A)  for  A  <  5 
was  not  performed  since  this  was  outside  of  the  experimental 
conditions;  the  computed  minimum  A  value  was  about  6. 
Figure  5.34  shows  that  the  least-squares-fit  curve  has  been 
extrapolated  for  A  >  18 .  Even  though  such  overextrapola- 
tions  are  not  generally  recommended,  it  appears  reasonable 
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for  this  situation.   The  justification  for  this  extapolation 
lies  with  the  diminishing  value  of  the  slope  m. 

Substituting  equation  (5.7)  into  (5.8)  yields: 


1  <j>2  (5.10) 

F2((j),A)  =  -[FxCO.AH  ~  m  T  ] 


In  order  to  modify  the  Fujii  et  al.  model,  three  assumptions 
were  made:  1)  heat  transfer  through  the  wire  is  negligible 
compared  to  that  through  the  interwire  space  for  the 
unflooded  portion  of  the  tube;  2)  heat  transfer  through  the 
flooded  portion  is  by  one-dimensionalconduction,  and  3)  in 
the  flooded  portion  of  the  tube,  heat  transfer  through  the 
condensate  between  the  wires  is  negligible  compared  to  the 
heat  flow  across  the  wires.  Based  on  assumption  (1)  above, 
the  heat-transfer  coefficient  for  the  unflooded  portion  can 
be  expressed  as  follows: 

.   hu       s      F2(p,A)  (5.11) 

hs    (s  +  Dw)   F2(0,0) 


where  Fo(0,0)  is  obtained  numerically  using  equation  (5.8). 
and  where  h„  is  the  heat- transfer  coefficient  for  the  smooth 
tube  predicted  by  Nusselt  theory  [27].  Also,  the  heat- 
transfer  coefficient  across  the  wires  in  the  flooded  portion 
of  the  tube  can  be  computed  by  equation  (5.12)  (derived  in 
Appendix  A) : 


1    rD        dx  (5.12) 

hf  -  d   /  W  2yU)  ,  Uw  -  -WT 
w  °    "T7-  + k 
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Now,  combining  equations  (5.11)  and  equation  (5.12),  equa- 
tion (5.13)  can  be  obtained  in  order  to  express  the  average 
heat- transfer  coefficient  for  the  entire  wire-wrapped  tube. 

h  -ha  (1  -  ~)  +hf  -  (7h^T  <5-13) 

Then  the  enhancement  ratio  is  given  by  equation  (5.14): 

E  -  "  (i  -  :)+rI:  _«  (5.14) 

0   hs       tt     hs  n   (s+Dw)  v 


In  order  to  fit  the  experimental  data,  a  value  of  0.02  was 
selected  for  the  constant  C  in  equation  (2.61).  Figure  5.35 
shows  the  experimental  data  and  the  values  calculated  from 
the  modified  Fujji  et  al.  model.  As  can  be  seen  from  this 
figure,  good  agreement  of  the  experimental  data  and  the 
predicted  values  exists.  However,  a  clear  trend  does  not 
exist  for  the  effect  of  wire  diameter  on  the  heat-transfer 
performance.  The  assumptions  made  during  this  study,  and 
other  possible  mechanisms  not  taken  into  consideration  owing 
to  the  complex  nature  of  this  problem,  such  as  convective 
effects  may  be  responsible  for  the  observed  comparison. 

H.   HEAT- TRANSFER  COEFFICIENT  FOR  COMMERCIAL  TUBES 

As  mentioned  in  Chapter  IV,  two  commercially  available 
finned  tubes  were  tested.  These  tubes  were  manufactured  by 
High  Performance  Tube,  Inc.,  and  had  fin  densities  of  1.1 
and  1.18  fins/mm  (tubes  51  and  49).  Figures  5.36  and  5.37 
show  the  variation  of  the  heat-transfer  coefficient  for 
these  tubes  with  heat  flux  under  vacuum  and  at  atmospheric 
pressure,  respectively.  Data  for  a  smooth  tube  (tube  50), 
are  also  shown.  The  best  performance  was  given  by  the  tube 
with   fin   density   of   1.1   fins/mm   under   both   pressure 
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conditions.  Visual  examination  of  the  condensation  process 
revealed  that  both  finned  tubes  were  fully  flooded.  This 
was  also  observed,  under  static  conditions.  Since  these 
tubes  had  a  high  density  of  fins  with  small  fin  height,  the 
surface  tension  effect  is  not  important  for  the  flow  of  the 
condensate  along  the  fin  side.  However,  the  surface  tension 
retained  the  condensate  between  the  fins,  so  these  tubes 
were  fully  flooded.  Therefore,  as  Table  III  shows,  although 
they  have  the  largest  area  ratio,  the  heat- transfer 
performance  is  worse  than  most  of  the  other  copper  finned 
tubes  tested  during  this  thesis  effort.  Nevertheless,  even 
though  these  tubes  were  fully  flooded,  the  normalized 
enhancement  ratio  (Eo/Ar)  was  greater  than  unity  in  both 
cases  (especially  for  atmospheric  pressure),  indicating  an 
enhancement  greater  than  first  due  to  an  area  increase. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

1.  The  use  of  fins  lead  to  significant  enhancement  of  the 
steam-side  heat-transfer  coefficient,  which  is  greater 
than  the  area  ratio  (finned  tube  area  to  smooth  tube 
area)  despite  condensate  retention  between  fins. 

2.  Enhancement  ratios  as  high  as  1.84  and  2.6  under 
vacuum  and  at  atmospheric  conditions,  respectively, 
were  realized  for  the  wire-wrapped  tube  with  a  wire 
diameter  of  0.5  mm  and  a  pitch  of  5.1  mm.  This  is  due 
to  the  surf ace- tension .effect  resulting  in  thinning  of 
the  condensate  between  wires. 

3.  The  Webb  et  al.  [31]  model  was  succesfully  modified  to 
predict  the  condensate  retention  angle  for  the 
wire-wrapped  tubes.  Maximum  error  of  15  percent  was 
found. 

4.  The  Fujii  et  al.  [32]  model  was  modified  to  predict 
the  vapor-side  "  heat-transfer  coefficient  for 
high-surface-tension  fluids  (i.e.,  water).  A  favorable 
agreement  between  the  modified  Fujii  et  al.  model  and 
the  experimental  data  was  found. 

5.  The  tube  with  a  "parabolic"  fin  profile  outperformed 
the  tubes  with  triangular,  trapezoidal  and  rectangular 
fin  shapes.  For  the  tube  with  "parabolic"  fins, 
enhancement  ratios  of  4.1  and  6.2  under  vacuum  and  at 
atmospheric  pressure,  respectively,  were  obtained. 

6.  For  tubes  with  spiral  triangular  fins,  the  optimum 
pitch  was  found  to  be  about  of  1.6  mm.  Enhancement 
ratios  of  3.9  and  6.1  under  vacuum  and  at  atmospheric 
pressure,  respectively,  were  obtained  for  this  tube. 

7 .  The  High  Performance  tubes ,  although  they  have  the 
highest  area  ratio  among  all  the  tubes  tested,  show 
poorer  performance  than  most  of  the  other  tubes.  This 
is  mainly  owing  to  the  high  condensate  retention 
angle.  Enhancement  ratios  as  high  as  3.5  and  5.2 
under  vacuum  and  at  atmospheric  pressure, 
respectively,  were  obtained. 

8.  The   enhancement   ratios   for   finned   tubes   are 

?roportional  to  the  tube  thermal  conductivity.  For 
he  tubes  with  different  thermal  conductivity,  the 
highest  enhancement  was  obtained  for  the  copper  tube, 
while  the  stainless  steel  tube  had  a  performance  even 
less  than  the  smooth  tube  under  vacuum  conditions. 
This  is  due  to  the  high  wall  thermal  resistance  of  the 
stainless  steel  tube. 
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B .   RECOMMENDATIONS 

1.  Attach  drainage  strips  on  some  of  the  existing  tubes 
to  investigate  the  effect  of  the  strips  on  the 
condensate  retention  angle  and  the  vapor- side 
heat-transfer  coefficient  and  compare  with  the  already 
existing  data. 

2.  Take  data  on  a  tube  with  a  more  nearly   parabolic  fin 

frofile  and  compare  them  with  the  existing  data  of 
ins   with   triangular,   trapezoidal   and   rectangular 
profiles . 

3.  Take  data  with  different  tube  diameters  to  investigate 
its  effect  on  the  heat-transfer  performance. 

4.  Take  data  with  different  vapor  velocities  to 
investigate  the  effect  of  vapor  shear  on  the 
heat-transfer  performance. 

5.  Take  data  using  different  fluids  to  study  the 
dependence  of  the  heat-transfer  coefficient  on  the 
fluid  properties. 
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APPENDIX  A 

DERIVATION  OF  THE  HEAT-TRANSFER  COEFFICIENT  OF  FLOODED 

REGION  FOR  WIRE-WRAPPED  TUBES 

Assuming  one-dimensional  heat  conduction  for  a  differen- 
tial element  dx ,  the  differential  heat  transfer  is  given  by 
equation  (A. 1) : 


AT 

dQ  =  ntoo  *  D°  dx  (A-1) 


where 


2y(x)   Dw  _  2y(x) 
ER(x)--  -Y—  +  -^T- (A. 2) 

where  y(x)  is  the  vertical  distance  between  the  wire  surface 
and  the  horizontal  tube  surface  (see  Figure  A.l),  and  is 
given  by  equation  (A. 3): 


2     2  !/2 
y(x)  =  Rv  ~  (Rw   "  x  }  (A'3) 

x  is  the  axial  coordinate  along  the  tube.   Substituting 
eqation  (A. 2)  for  £R(x)  into  equation  (A.l)  and  integrating 


2  d   /2  ^x 

Q   =    tt   DQ    AT   D-        /w        2yCxT Uw  -   ggg  (A'4) 

kf  kw 


Also  heat  transfer  is  given  by  equation  (A. 5) 

141 


w 


Q  -  h  it  DQ  -x—  AT 


(A. 5) 


Combining  equations  (A. 4)  and  (A. 5),  the  heat-transfer  coef- 
ficient can  be  expressed  as: 


1    D         dx 

hf  -  d  i  w  "^y^n — uw  -  2KD 
w  °       ir  +  — i- 


(A. 6) 


Figure  A.l   Condensate  Film  Profile  for 
Fully  Flooded  Wire-Wrapped  Tube. 
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APPENDIX  B 
LISTING  OF  RAW  DATA 

The  following  pages  contain  raw  data  obtained  for  tubes 
number  6,  17,  36,  38  and  45  thru  71  under  vacuum  conditions 
and  at  atmospheric  pressure. 
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APPENDIX  C 
UNCERTAINTY  ANALYSIS 

Experimentally  determined  quantities  are  always  associ- 
ated with  uncertainties  owing  to  the  measuring  device  accu- 
racy, calibration  of  the  device,  and  the  operator's 
experience.  During  this  thesis  effort,  numerical  data  were 
taken  and,  together  with  theoretical  formulation,  the  steam- 
side  heat- transfer  coefficients  were  calculated.  Since  the 
devices  used  during  this  experiment  to  read  steam  tempera- 
ture, inlet  and  outlet  cooling  water  temperatures,  flowrate 
of  the  cooling  water  in  finally  computing  the  steam-side 
heat- transfer  coefficient,  the  final  result  may  be  distorted 
due  to  the  uncertainty  propagation  during  calculations.  In 
cases  where  the  final  results  show  large  uncertainties,  it 
may  be  unwise  to  accept  the  experimental  results.  The 
uncertainty  on  a  computation  can  be  determined  using  the 
following  equation  proposed  by  Kline  and  McClintok  [35] 
shown  below: 


Wo  - 


3R 


dX, 


w. 


1/2 


(C.l) 


where 


W^  is  the  uncertainty  of  the  desired  dependent 
variable 

xl»  x2 '  ••••»  xn  are  t^ie  measured  (independent) 
variables 

Wi,  Wo,  ....,  w   are  the  uncertainties  in  the  measured 
variables 
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Using  program  "UNA7"  which  is  listed  at  the  end  of  this 
appendix  the  uncertainties  associated  with  various  quanti- 
ties during  this  investigation  were  obtained.  Also,  listed 
in  this  appendix  are  some  of  the  selected  uncertainty  evalu- 
ations. A  complete  discussion  on  the  uncertainty  analysis 
used  for  this  experiment  is  given  by  Georgiadis  [5] . 
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FILE  NAME 
REVISED  : 


UNA7 

February  28.  1986 


COM  /Cc/  Cl7) 

DIM  E(4) 

DATA  0. 10086091 ,25727.94369.-767345.8295.78025595.81 

DATA  -9247486589, 6. 97688E+ 11 , -2.661 92E+1 3 .3.94078E+1 4 

READ  C<*> 

PRINT 

PRINTER  IS  701 

PRINT  USING  "1  OX." "DATA  FOR  THE  UNCERTAINTY  ANALYSIS:*"" 

PRINT 

BEEP 

INPUT  "ENTER  FILE  NAME".File$ 

PRINT  USING  "10X,""File  Name:  " "  .  1 2A"  ;F  i  leS 

BEEP 

INPUT  "ENTER  DATA  SET  NUMBER  FOR  UNCERTAINTY  ANALYSIS" , Ids 

BEEP 

INPUT  "ENTER  PRESSURE  CONDITION  (0=V,1=A)  ".Pre 

Prc=Prc+1 

BEEP 

INPUT  "ENTER  Ci'.Ci 

ASSIGN  SFile  TO  FileS 

ENTER  C*File:Ifg.Inn 

IF  Ifg-0  THEN  ENTER  *File:Dd 

IF  Ifg-1  THEN  ENTER  »Fi le:Dd.Dd, Dd 

FOR  1=1  TO  Ids 

ENTER  @File:Bvol .Bamp , Vtran.Etp,E( * ) .Fm.Tci . Tco .Phg.Puater 

NEXT  I 

Emf=E(0) 

IF  Prc=1  THEN 

BEEP 

PRINT  USING  "10X. ""Pressure  Condition:    Vacuum  (11  kPa)""" 

ELSE 

PRINT  USING  "10X, ""Pressure  Condition:    Atmospheric  (101  kPa) 

END  IF 

BEEP 

PRINTER  IS  1 

PRINT  USING  "4X, 

PRINT  USING  "6X 

PRINT  USING  "6X, 

PRINT  USING  "SX, 

INPUT  Itt 

IF  ltt=0  THEN 

BEEP 

INPUT  "SELECT  (0-THIN,  1 -THICK) ", Iut 

END  IF 

PRINTER  IS  701 

IF  Itt-0  THEN 

Di=.0127      !  Inside  diameter  of  test  tube 

Kc=385 

Dkc=10 

IF  Iut-0  THEN 

Do=.0137 

ELSE 

Do=. 01905     !  Outside  diameter  of  test  tub* 

END  IF 


"Select  Material  Code:""" 
"'0  Copper     1  Stainless  steel 

"2  Aluminum   3  90:10  Cu-Ni 

"4  HFTI""" 
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1310  END  IF 

1315  IF  Itt-1  THEM 

I320  Kc=16 

1325  Dkc=1 

1330  Di=. 01245 

1335  Do-. 0145 

1340  END  IF 

1345  IF  Itt=2  THEN 

1350  Kc=167 

1355  Dkc=5 

1360  Do=.0137 

1365  Di-.0127 

1370  END  IF 

1375  IF  Itt=3  THEN 

1380  Kc=45 

1385  Dlcc  =  2 

1390  Di-,0127 

1395  Do=.0137 

1400  END    IF 

1405  IF  Itt=4  THEN 

1410  Kc=38S 

1415  Dkc=10 

1420  Di=.0156 

1425  Do=.0175 

1430  END  IF 

1435  D1=. 01905 

1440  D2-. 01587 

1445  IF  Itt=4  THEN  D2=. 01905 

1450  PRINTER  IS  701 

1455  Ts=FNTvsv(Emf ) 

1460  PRINT  USING  " 1 0X. " "Steam  Temperature  =  "B.3D.2D,""        (Deg 

C) :Ts 

1465  PRINT  USING  " 1 0X. " "Water  Flow  Rate  (%)  =  "".3D.2D";Fm 

1470  Dtci=.01 

1475  Dtco-.OI 

1480  BEEP 

1485  Demt-I.OE-G 

1490  Dts  =  SQR(  (  l  C  C  1  )+2+C(2)*Emf  +  3*C(3)-*Emf  2+4*C( 4 )*Emt  *3>*Demf  )  "2) 

1495  T=(Tci+Tco)/2   !  FILM  TEMPERATURE 

1500!  UNCERTAINTY  IN  THE  COOLING  WATER 

1505  Drho=.S  !  ERROR  IN  WATER  DENSITY 

1510  Dmt=.0Q44  !  ERROR  IN  MASS  FLOW  RATE 

1515  Rho=FNRho(T)  !  WATER  DENSITY 

1520  Mf=1 .04805E-2+6.80932E-3*Fm   !  MASS  FLOW  RATE  OF  COOLING  WATER 

1525!  CORRECT  MF  FOR  THE  TEMPERATURE  EFFECT 

1530  Mf=Mf+(1 .0365-1 . 96644E-3*Tci+S.2S2E-6*Tci *2> /. 995434 

1535  Ai=(PI*Di'2)/4  !  TUBE  INSIDE  CROSS  SECTION  AREA 

1540  Ddi-. 000025 

1545  Dai=PI*Di*Ddi/2  !  ERROR  OF  INSIDE  TUBE  CROSS  AREA 

1550!  COMPUTE  THE  WATER  VELOCITY 

1555  Vu=Mf/(Rho*Ai)  !  WATER  VELOCITY 

1560  PRINT  USING  "10X, ""Water  Velocity  =    "".Z.DD,""        Cm/ 

3)""":Vu 

1565!  CORRECT  OUTLET  WATER  TEMP.  FOR  THE  MIXING  CHAMBER  EFFECT 

1570  IF  Inn-1  OR  Inn=5  THEN  Tco=Tco- . 004*Vu - 2 

1575  IF  lnn=0  THEN  Tco-Tco- < - . 00138+. 00T*Vu*2> 

1580  T=(Tci+Tco)*.S  !  FILM  TEMPERATURE 

1585!  COMPUTE  THE  ERROR  IN  WATER  VELOCITY 

1590  Dvw=Vw+SQR( (Dmf/Mf ) " 2+ ( Drho/Rho ) ' 2+( Dai/Ai ) "2 > 

1595!  UNCERTAINTY  IN  THE  REYNOLDS  NUMBER 
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1600  Mu=FNMw(T)  !  WATER  VISCOSITY 

1G0S  Dmu=6.E-6  !  ERROR  OF  WATER  VISCOSITY 

1610  Re=(Rho*Vu*Di )/Mu 

161 S  Dre=Re*SOR< (Drho/Rho) *2+(Dvu/Vu) '2+(Ddi/Di > "2+ ( Dmu/Mw) '2) 

1620!  UNCERTAINTY  IN  THE  HEAT  TRANSFERRED 

1625  Cpw=FNCpw(T) 

1630  Q=Mf*(Tco-Tci)*Cpu 

163S  Dcpu»=8 

1640  Dq  =  Q*SQR((Dmf/Mf ) -2+( (Dtco/( Tco-Tci ) )) *2+( (Dtci/( Tco-Tci ) ) )  "2+ ( Dcpu/Cpu) "2 

) 

1645!  UNCERTAINTY  IN  THE  HEAT  FLUX 

16S0  Dl= .0005  !  ERROR  IN  TUBE  LENGHT 

16SS  Ddo=.00002S 

1660  L=.  13335  !  CONDENSING  TUBE  LENGTH 

1665  Qp=Q/(PI*Do*L)  !  HEAT  FLUX 

1670  PRINT  USING  "10X,"MHeat  Flux  ■    "",Z.3DE."U   (W/m'2)" 

"":QP 

1675  PRINT  USING  " 1 0X ." "Tube-metal  thermal  conduc.     =  "".3D.D,""        (W/m. 

K) ;Kc 

1680  PRINT  USING  " 1 0X. ""Sieder-Tate  constant  -    "".Z.4D":Ci 

1685  Dqp=Qp*SQR( (Dq/Q) *2+( Ddo/Do ) *2+(Dl/L) *2) 

1690  Lmtd=(Tco-Tci)/LOG( (Ts-Tci )/< Ts-Tco) ) 

1695  Uo=Qp/Lmtd  !  OVERALL  HEAT  TRANSFER  COEF. 

1700  A1«Dts*(Tci-Tco)/((Ts-Tci)*< Ts-Tco )*L0G< ( Ts-Tci ) /( Ts-Tco)  )  ) 

1705  A2  =  Dtci/( (Ts-Tci)*L0G( ( Ts-Tci )/( Ts-Tco)  )  ) 

1710  A3  =  Dtco/( ( Ts-Tco )*L0G( (Ts-Tci )/( Ts-Tco) )  ) 

1715  Dlmtd=Lmtd*S0R(A1 * 2+A2*2+A3" 2 ) 

1720  Duo=Uo*SQR( (Dqp/Qp) "2+(Dlmtd/Lmtd) '2) 

1725  M=Mu 

1730  T1=(T+273.15)/273.15 

1735  Ku=FNKw(T1) 

1740  Ac=0.    !  INTERSCEPT  FROM  SIEDER  PROGRAM 

1745  L1=. 060325  !  LENGTH  OF  UNFINNED  LEFT  PART  OF  TUBE 

1750  L2=. 034925  !  LENGTH  OF  UNFINNED  RIGHT  PART  OF  TUBE 

1755  Pr=Cpw*Mu/Kuj 

1760  Muu=FNMutj(T) 

1765!  UNCERTAINTY  OF  INSIDE  HEAT-TRANSFER  COEFF. 

1770  Cf-1. 

1775  Hi=(Ku/Di)*(Ci*Re*.8*Pr" .333*Cf+Ac) 

1780  Dti=Q/(PI*Di*(L+L1*Fe1+L2*-Fe2)*Hi) 

1785  Cfc=(Muu/FNMuw(T  +  Dti) ) *  .14 

1790  IF  ABS<<Cfc-Cf )/Cfc)>.01  THEN 

1795  Cf=(Cf+Cfc)*.S 

1800  GOTO  1775 

1805  END  IF 

1810  PI =PI*( Di+DI  ) 

1815  B1-(D1-Di)*PI*(Di+D1)*.5 

1820  M1=(Hi*P1/(Kc*B1 ) )' .5 

1825  P2=PI*(Di+D2) 

1830  B2=(D2-Di)*PI*(Di+D2)*.5 

1835  M2=(Hi*P2/(Kc«B2> )  '.5 

1840  Fe1=FNTanh(M1*L1 )/<M1*L1 ) 

1845  Fe2=FNTanh<M2*L2)/(M2+L2) 

1850  Dtc=Q/(PI*Di*(L+L1*Fe1+L2*Fe2)*Hi ) 

1855  IF  ABS((Dtc-Dti)/Dtc)>.01  THEN  1775 

1860  Dlcw-.OOIO  !  ERROR  IN  WATER  THERMAL  CONDUCTIVITY 

1865  Dci=.0005  !  ERROR  IN  SIEDER-TATE  COEFFICIENT 

1870  Dpr=.0S  !  ERROR  IN  PRANDTL  NUMBER 

1875  Dcf=8.E-6 

1880  A4=.14*Dcf/Cf 
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1885 

A4) 

1890! 

1895 

1900 

1905 

i)  ))' 

1910 

1915 

1920 

1925 

1930 

1935 

1940! 

1945 

1950! 

1955 

I960! 

1965 

1970! 

1975 

1980! 

1985 

1990! 

1995 

2000! 

2005 

2010! 

2015 

2020 

2025 

2030 

2035 

2040 

2045 

2050 

2055 

2060 

2065 

2070 

2075 

2080 

2085 

2090 

2095 

2100 

2105 

2110 

2115 

2120 

2125 

2130 

2135 

2140 

2145 

2150 

2155 

2160 

2165 

2170 


Dhi=Hi«SQR< (Dku/Kw) *2+( Ddi/Di ) "2+( .8*0  re/Re) "2  +  ( . 333*Dpr/Pr ) "2+ < Dci/Ci ) '2  + 


UNCERTAINTY  OF  OUTSIDE  HEAT-TRANSFER  COEFF. 
Ru=Do*L0G(Do/Di)/(2*Kc)  !  WALL  RESISTANCE 
Ho=1/( ( 1/Uo)-(Do*L/(Di*(L+L1+Fe1+L2*Fe2)*Hi) )-Rw) 
Dru  =  Rw*SQR( (Ddo/Do) *2+( Dkc/Kc ) ' 2+ ( Ddo/ ( Do*L0G( Do/D  i  )  )  ) 
2) 
A5=1/Uo-Ru-<Do+L/(Di*(L+L1*Fe;+L2*Fe2)*Hi) ) 
A6=Duo/(Uo'2*A5) 
A7=Dru/AS 

A8=( (Do/(Di*Hi) )*(Dhi/Hi ) ) /AS 
PRINT 

Dho=Ho*SQR(A6"2+A7'2+A8'2) 
CALCULATE  THE  Z  UNCERTAINTY 


2+(Ddi/(Di*L0G(Do/D 


IN 

Prho=Dho*100/Ho 

CALCULATE  THE  7.  UNCERTAINTY  IN 
Prre=Dre*100/Re 

CALCULATE  THE  7.  UNCERTAINTY  IN 
Prmf=Dmt*100/Mf 

CALCULATE  THE  7.  UNCERTAINTY  IN 
Prqp=Dqp+100/Qp 

CALCULATE  THE  7.  UNCERTAINTY  IN 
Prlmtd  =  Dlmtd*100/Lintd 
CALCULATE  THE  7.  UNCERTAINTY  IN 
Prrw  =  Dru+100/Rij 

CALCULATE  THE  7.  UNCERTAINTY  IN 
Pruo=Duo+100/Uo 

CALCULATE  THE  7.  UNCERTAINTY  IN 
Prhi=Dhi*100/Hi 
PRINT 

10X, 


Ho 

REYNOLDS  NUMBER 

MASS  FLOW  RATE 

HEAT  TRANSFER 

LMTD 

Ru 

OVERALL  HEAT  TRANSFER  COEF. 

INSIDE  HEAT  TRANSFER  COEFF. 


USING 

USING 

USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 


10X. 
10X. 

1  ox. 

10X. 

1  ox. 

10X. 
10X. 
10X, 
10X. 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

END 

DEF  FNMuu(T) 

A=247.8/(T+133.15) 

Muu=2.4E-5*10"A 

RETURN  Muu 

FNEND 

DEF  FNTanh(X) 

P=EXP(X) 

Q-EXP(-X) 

Tanh=(P-Q)/(P+Q) 

RETURN  Tanh 

FNEND 

DEF  FNKu(T1 ) 

Ku=-.92247+T1*(2.839S 

RETURN  Ku 

FNEND 

DEF  FNMu(T) 

A=247.3/(T+133.15) 


UNCERTAINTY  ANALYSIS: 


VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  Number.  Re 
Heat  Flux,  q 
Log-Mean-Tem  Diff.  LMTD 
Wall  Resistance,  Ru 


PERCENT  UNCERTAINTY" 


iPrmf 


Z.2D. 
Z.2D.":Prre 
Z.2D.";Prqp 
,Z.2D":Prlmtd 


Overall  H.T.C. 


Water-Side 
Steam-Side 


Uo 


Hi 
Ho 


".DD.2D. 
"  .DD.2D, 
.3D.2D.' 

.3D. 2D,' 


";Prrw 
":Pruo 
;Prhi 
;Prho 


■T1*( 1 . 8007-T1 *( . 52577-. 07344*T1 ))) 
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2175  Mu-2.4E-5*10"A 

2180  RETURN  flu 

2135  FNEND 

2190  DEF  FNRho(T) 

2195  Rho=999.S2946+T*( . 01269-T*(5. 482513E-3-T*1 .234147E-5)  ) 

2200  RETURN  Rho 

2205  FNEND 

2210  DEF  FNCpw(T) 

2215  Cpu=(4.21 1208S8-T*(2.26826E-3-T*(4.42361E-S+2.71428E-7*T) ) )+l 000 

2220  RETURN  Cpu 

2225  FNEND 

2230  DEF  FNTvsv(Emf) 

2235  COM  /Cc/  C(7) 

2240  T-C(O) 

2245  FOR  1=1  TO  7 

2250  T-T+C(I)*Emf"  I 

2255  NEXT  I 

2260  RETURN  T 

2265  FNEND 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name : 

Pressure  Condi t ion: 

Steam  Temperature 

Hater  Flow  Rate  i%) 

Water  Velocity 

Heat  Flux 

Tube-metal  thermal  conduc. 

Sieder-Tate  constant 


F06V145 

Vacuum  ( ! 1 


kfa  > 
46.55 
SO.  00 

4 .  39 

3.329E+05 


38c 


3  .  U 

0.0658 


! Deq  C ) 

i.m/'s  ) 
(W/m*2) 
(W/m.K) 


UNCERTAINTY  ANALYSIS: 

VARIABLE 

Mass  Flow  Rate.  Md 
Reynolds  Number.  Re 
Heat  Flux ,  q 
Log-Mean-Tem  Ditt,  LMTD 
Wall  Resistance,  Ru 
Overall  H.T.C. .  Uo 
Water-Side  H.T.C. ,  Hi 
Steam-Side  H.T.C,  Ho 


PERCENT  UNCERTAINTY 


0 

,79 

1 

1 

1  0 
1  c. 

1 
1 

,53 

1 

24 

67 

]       ( 

97 

1 

1  ■_  w> 

w>     1 

10 

DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


F  i  le  Name:  F06V1 45 

Pressure  Cond 1 1  ion :    Vacuum 

Steam  Temperature 

Water  Flow  Rate  (7) 

Water  Ve loci ty 

Heat  Flux 

Tube-metal  thermal  conduc. 

Sieder-Tate  constant 


kra) 
48.  39 
20.00 

1.16 

1 .842E+0S 
385.0 

0.0658 


(Deg  C) 

(m/s) 

(W/m'2) 
(W/m.K) 


UNCERTAINTY  ANALYSIS: 

VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux,  q 
Log-Mean-Tem  Ditf.  LMTD 
Wall  Resistance,  Rw 
Overall  H.T.C. ,  Uo 
Water-Side  H.T.C. .  Hi 
Steam-Side  H.T.C,  Ho 


PERCENT  UNCERTAINTY 


3 . 0 1 


*  J 

.10 

6 

.60 

-~> 

.67 

-j 

.  15 

~> 

.63 

>o 

.05 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name: 

Pressure  Condition: 

Steam  Temperature 

Water  Flow  Rate  (%) 

Na ter  Velocity 

Heat  Flux 

Tube-metal  thermal  conduc. 

Sieder-Tate  constant 


F06A226 

Atmospher ic  (101 
=   99.96 


Pa) 


80 

4 

1 

38S 

0 


00 

40 

038E+06 
0 

0630 


(Deg  C) 

(m/s) 
(W/m*2) 

(N/m.K) 


UNCERTAINTY  ANALYSIS: 

VARIABLE 

Mass  Flou  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux ,  q 
Log-Mean-Tem  Diff,  LMTD 
Wall  Resistance,  Ru 
Overall  H.T.C. .  Uo 
Water-Side  H.T.C. ,  Hi 
S team-Side  H.T.C,  Ho 


PERCENT  UNCERTAINTY 

0.79 
1.12 

0.99 
0.40 
2.67 
1  .07 
1  .25 
3.22 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name : 
Pressure  Condition: 
Steam  Temperature 


F06A226 
Atmospheric  ( 
=   99 


Water  Flou  Rate 

Water  Veloci ty 

Heat  Flux 

Tube-metal  thermal  conduc. 

Sieder-Tate  constant 


20 
1 

6 
385 

0 


01  kPa) 

99 

00 

16 

036E+05 

0 

0630 


(Deg  O 

(m/s) 

(W/m'2) 
(W/m.K) 


UNCERTAINTY  ANALYSIS: 

VARIABLE 

Mass  Flou  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux,  q 
Log-Mean-Tem  Diff,  LMTD 
Wall  Resistance,  Ru 
Overall  H.T.C. .  Uo 
Water-Side  H.T.C. ,  Hi 
Steam-Side  H.T.C. ,  Ho 


PERCENT  UNCERTAINTY 

3.00 
3.1  1 
3.04 
0 .  1 8 
2.67 
3.04 
2.64 
19.36 
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